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Chapter #tie 

INTRODUCTION 
FLUID PROPERTIES 


1.1. Definitions : 

A fluid is defined as a substance that deforms continuously whilst acted upon by any force 
tangential to the area on which it acts. Such a force is termed a shear force, and the ratio of the 
shear force to the area on which it acts is known as the shear stress. 

Fluid Mechanics: Fluid mechanics may be defined as that branch of engineering science that 

deals with the behavior of fluid under the condition of rest and motion 

Fluid mechanics may be divided into three parts: Statics, Kinematics, and Dynamics 

. Statics Deals with fluid at rest in equilibrium state, no force no acceleration. 

. Kinematics Deals With flow behaviors of fluid like velocity, acceleration and flow patterns. 

. Dynamics Deals with the effects of flow behaviors on fluid surroundings like forces and 
momentum exchange. 


Contact area. Shear stress 

A \ / t-F/A Force, F 



Shear 
strain, a 
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1.2. Fluid Properties : 

1. Mass Density (p) (Specific Mass) : the ratio of the mass of a given amount of the 
substance to the volume it occupies. For liquids, this definition is generally satisfactory 


= — (kg/m 3 ).... (1.1) 


where m : mass of the fluid (kg) , v : volume (m 3 ) 


Mean Density : the ratio of the mass of a given amount of a substance to the volume that this 
amount occupies. If the mean density in all parts of a substance is the same then the density is 
said to be uniform. 

Relative Density: is the ratio of the density of a substance to some standard Relative density 
density. 

2. Weight Density (y) (Specific Weight) : is defined as the weight per unit volume at the 
standard temperature and pressure, 


= -= p*g (N/m 3 ).... (1.2) 


where w : weight of the fluid (kg) 


, v : volume (m 3 ) 


g : gravitational acceleration = 9.81 m/s 2 


y of water = 9810 N/m 3 at 4°C and 1 Bar 


3. Specific Volume (v): It is defined as a volume per unit mass of fluid, 



Its unit are m 3 /Kg. 


4. Specific Gravity (Sg): It is defined as the ratio of the specific weight of the fluid to the 
specific weight of a standard fluid 
For liquids the standard fluid is pure water at 4°C, 



For Gasses the standard fluid is air 
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example ft) : Calculate the Specific weight, specific mass, specific volume and specific gravity 
of a liquid having a volume of 6m 3 and weight of 44 kN. 

Solution: 


W=44 kN 
W 



, V= 6 m3 
44 

— = 7.333 kN/m 3 
6 

7.333 x 1000 

— — = 747.5 kg/m 3 

9.81 


1 _ 1 
p ~~ 747.5 


0.00134 m 3 /kg 


Sg = 


Y liquid 

Ywater 


7.333 

9.81 


0.747 


example (2): Calculate the Specific weight, specific mass, specific volume and density of chlorine 
gas(cl2) at 80° F and pressure of 100 psia(Ibs per square inch, absolute) 

Solution: 

Atom. Wt. chlorine = 35.45 mol. Wt. = 2*35.45 = 70.9 
R = 49709/70.9 = 701.1 ft . lb / slug °R 
y = 32.2 * 100 * 144 / 701.1 * 540 = 1.225 lb/ft 3 
specific volume = 1/1.225 = 0.815 ft 3 /lb 
p = 1.225/32.2 = 0.0381 slug/ft 3 


Example (3) : Calculate the Specific weight, specific volume and density of carbon dioxide 
gas(C02) at 100°C and atmospheric pressure of 101.3 kN/m 2 (kilo Newton per square meter, 
absolute) 

Solution: 

Atom. Wt. carbon = 12 , Atom. Wt. oxygen =16 mol. Wt. = 2*16+12 = 44 
R = 8313/44 = 189 ft . lb / slug °R 
y = 9.81 * 101.3 * 10 3 / 189 * 373 = 1.225 lb/ft 3 
specific volume = 1/14.1 = 0.071 m 3 /N 
p = y/g = 14.1/9.81 = 1.44 kg/m 3 
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5. Elasticity ( Compressibility) : 

All matter is to some extent compressible. That is to say, a change in the pressure applied 
to a certain amount of a substance always produces some change in its volume. Although the 
compressibility of different substances varies widely, the proportionate change in volume of 
a particular material that does not change its phase (e.g. from liquid to solid) during the 
compression is directly related to the change in the pressure. 

Bulk modulus of The degree of compressibility of a substance is characterized by the modulus 
of elasticity, E, which is defined by the equation: 


There is small pressure disturbance, that travels as a wave of increased (or decreased) density 
and pressure, moves at a celerity, a , given by : 



F 



Strain V/Vi 


1 



Where a is frequently called the sonic or acoustic velocity because sound, a small pressure 
disturbance, travels at this velocity. 
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example (4) : Air at 60° F and 14.7 psia is compressed isentropically so that its volume is reduced 
50%. Calculate the final pressure and temperature and the sonic velocities before and after 
compression. 

Solution: 

yi = 32.2 * 14.7 * 144 / 1715 * 520 = 0.0763 lb/ft 3 

y 2 =2* 0.0763 = 0.1526 lb/ft 3 

p 2 * 144 14.7 * 144 

(0.1526) 1 - 4 “ (0.0763) 1 - 4 


P 2 = 38.8 psia 


32.2 * 38.8 * 144 
1715 *T 2 


0.1526 


T 2 = 687° R (227° F) 


ai = Vl.4 * 1715 * 520 = 49V520 = 1118 fps 
a 2 = 49 V687" = 1285 fps 


6. Viscosity : 

The most important of fluid properties is viscosity, which relates the local stresses in a moving 
fluid to the strain rate of the fluid element. When a fluid is sheared, it begins to move at a strain 
rate inversely proportional to a property called its coefficient of viscosity. Consider a fluid element 
sheared in one plane by a single shear stress, as in figure below. The shear strain angle will 
continuously grow with time as long as the stress is maintained, the upper surface moving at 
speed u larger than the lower. Such common fluids as water, oil, and air show a linear relation 
between applied shear and resulting strain rate 



from the geometry of figure below, we see that: 


tan d0 = — ....(1.8) 

dy 


5 
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In the limit of infinitesimal changes, this becomes a relation between shear strain rate and 
velocity gradient : 


d6 dv 

1 • • • • 

at ay 



From equation (1.7), then, the applied shear is also proportional to the velocity gradient for the 
common linear fluids. The constant of proportionality is the viscosity coefficient , p : 




Therefore, the coefficient of dynamic viscosity: 

(N. s/m 2 or Pa .s or 



kg/m .s) 


Kinematic Viscosity , v = the ratio of dynamic viscosity to mass density 


1 


The unit Stoke (St) is also used where 


1 st = 10' 4 m 2 /s (1 st = cm 2 /s) 


6 


I 
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• The fluid is non-Newtonian if the relation between shear stress and shear strain rate is 

non-linear 


• Typically, as temperature increases, the viscosity will decrease for a liquid, but will 
increase for a gas. 


example (5) : Oil has dynamic viscosity (ju = 1.0 * 10-3 Pa.s) filled the space between two 
concentric cylinders, where the inner one is movable and the outer is fixed. If the inner and 
outer cylinders has diameters 150mm and 156mm respectively and the height of both cylinders 
is 250mm, determine the value of the torque (T) that necessary to rotate the internal cylinder 
with 12 rpm? 

Solution: 

v = IE1H x 2ttx = — x In x 0.075 = 0.09425 m/s 
60 60 

T = xChr 

r = T ( 2 7tt xhx r)~ l = T ( lahr 2 ) _1 

dv Av 

T = i J ~r = i l ir 

dy A y 

, 0.09425 n H r m N 

r = 10 3 x = 31.41667 x 10 2 — 

0.003 m 2 

T = 31.41667 x 10“ 2 x 2 x .075 x n x .25 x .075 = 2.7 x 10“ 4 N. m 
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€xamole (6) : Oil has a density of 580 kg/m 3 flow through a pipe its diameter 200 mm. If it is 
known from the pressure calculations for a certain length of the pipe that the shear stress at 
the pipe wall equal to 0.07 N/m 2 , and its known from the velocity calculations through a certain 
cross section of the pipe that the velocity profile equation is : 

V = 1-lOOr 2 


Where the velocity dimension is in m/s and the distance from the centre of pipe r in m. If the 
flow is laminar, calculate the kinematic viscosity for the oil? 

Solution: 

For laminar flow : 

dv 

T = LI 

dy 

the distance from the centre of pipe is 
r = R-y 

Where R is the radius of the pipe, y is the distance from the pipe wall toward the pipe centre. 

Where dr = -dy and then the above equation become: 

dv 

T = -jU — 


And from the section of velocity distribution, the strain in any point equal to : 

— = -200r 
dr 

then the shear stress in any distance from the pipe centre is expressed as : 
t = 200 pr 

whereas z at the wall equal to 0.07 N/m 2 then by substitution this value in the last equation we obtain 
follows: 


0.07 N.s 


0.0035 Pa.s 


8 


and the kinematic viscosity is : 
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y = ^ = 0.0035^.m- i =41176xl0 ^A^ = 41176xl0 _ 8 m^ 


p 850kg.m 


-3 


kg 


example (7) : The velocity distribution for flow over a plate is given by u= 2y + y 2 where u is the 
velocity in m/s at a distance y meters above the plate surface. Determine the velocity gradient 
and shear stresses at the boundary and 1.5m from it. Take dynamic viscosity of fluid as 0.9 
N.s/m 2 

Solution: 

Given u= 2y + y 2 ^ = 2 — 2y 


1 ) Velocity gradient, — 
At boundary :at y=0 
At y =0.1 5m 


du „ _i 

— =2 s 1 

dy 

-=2-2x0.15 = 1.7s -1 

dy 


2) Shear stresses, x: 

du N 

T)y= o = p— = 0.9 x 2 = 1.8 — 
y dy m z 

du 

T )y=o.i 5 = jU— = 0.9 x 1.7 = 1.53 N /m 2 


Example (8) : A cylinder 0.28 ft. in radius and 2 ft. in length rotates coaxially inside a fixed 
cylinder of the same length and 0.3 ft radius. Glycerin (p = 0.031 lb.sec/ft 2 )fills the space 
between the cylinders. A torque of 0.5 ft.lb is applied to the inner cylinder. After constant 
velocity is attained, calculate the velocity gradient at the cylinder walls, the resulting rpm, and 
the power dissipated by fluid resistance. Ignore end effects. 


Solution: 




9 


I 
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Refer to figure 2 above for a definition sketch. The torque T is transmitted from inner cylinder to outer 
cylinder through the fluid layers; therefore ( r is the radial distance to any fluid layer and / is the length 
of the cylinders) 


T = -x(2 7tr*/)r ; t = — 
Consequently : 

d G) T 

dr 
and 


2nr 2 l 


M r d (-) 

dr 


2n filr 3 


( v )~ ~ T 

rr Q dr 

1 

^ 1 

1 

l 

\rj 2nfx l ' 

J ri r 3 

ri 4n /xl 

. ' o 

" ?. 


V = °' 28 * 0 ' 5 f = 0.296 ft/sec 

471*0.031*2 L (0.3) 2 (0.28) 2 J 

The power dissipated in fluid friction is : 

P = 2 7i n Ti V/ = Too 


Where to = V / n = 1.06 sec 1 is the rotational speed (radians/sec) 

Numerically : 

= ~ T + — = —16.4 + 1.1 = —15.3 ft / sec / ft 

2nn l rf ri 

= ~ T , + —= -14.2 + 0 = -14.2 ft /sec /ft 

2nfx l r£ r Q 
(^) * 60 = 10.1 
P = 0.5 * 1.06 = 0.53 ft.lb/sec = 0.00096 hp 

This power will appear as heat, tending to raise the fluid temperature and decrease its viscosity; 
evidently a suitable heat exchanger would be needed to preserve the steady state conditions given. 

When h = r 0 - n — H) , then h / r, is small and dv/dr — W / h 

As a consequence, v/r is much less than dv/d] * : 


^dr Jl 


A 

Mr 0 


rpm = 
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dv _ n V 


Assuming this linear profile case for an approximate calculation gives : 


h = 0.02 ; 


h 


0.07 


T 

V = — 


27TjU l 



1.28 


0.0784 


0.02 


= 0.327 ft / sec. = 11.2 rpm 


Because these results differ from the former by almost 11%, the approximation is not satisfactory 
in this case. 


7. Surface Tension ( Capillarity) : 

The apparent tension effects that occur on the surfaces of liquids, when the surfaces are in 
contact with another fluids or a solid, depend fundamentally upon the relative sizes of 
intermolecular cohesive and adhesive forces. Although such forces are negligible in many 
engineering problems, they may be predominant in some such as: 

1. The capillary rise of liquids in narrow spaces. 

2. The mechanics of bubble formation. 

3. The breakup of liquid jets. 


A molecule 
on the surface 




-A molecule 
inside the 
liquid 


Surface tension is a properties of certain fluid-fluid interface 


water-air 0.075 N/m at 20°C Water-air .... 0.056 N/m at 100°C 


mercury- air ... 0.1 N/m 


li 
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nK 2 A p 


2nRY 



pressure inside water droplet: 

let P= The pressure inside the drop 
d= Diameter of droplet 

g= Surface tension of the liquid (water-air interface) 
From sectional free body diagram of water droplet we have 



1 . 

2 . 

3 . 


AP between inside and outside = P-0 =P 
Pressure force =P X - d 2 , and 

4 

Surface tension force acting around the circumference= <r xnd, 
under equilibrium condition these two forces will be equal and opposite, i.e. 
P x-d 2 = o xnd 

4 


axnd 

4a 

D 

Jd 2 

4 

= — ....(1.13) 

a 


From this equation we show that (with an increase in size of droplet the pressure intensity 
is decreases) 


12 
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Notes : 

Property of exerting forces on fluids by fine tubes and porous media, due to both 
cohesion and adhesion 

Cohesion < adhesion, liquid wets solid, rises at point of contact 
Cohesion > adhesion, liquid surface depresses at point of contact, non-wetting fluid 
The contact angle is defined as the angle between the liquid and solid surface 
Meniscus: curved liquid surface that develops in a tube 


weight of fluid column = surface tension pulling force 


a(nR 2 h) = 2nRo cos0 ....(1.14) 


» 2 (TCOS0 

h = .... (1.15) 

pgR 


N 


Liquid 


ju 


2R 


IV 


lirRa, 




• Expression above calculates the approximate capillary rise in a small tube 

• The meniscus lifts a small amount of liquid near the tube walls, as r increases this amount 
may become insignificant 

• Thus, the equation developed overestimates the amount of capillary rise or depression, 
particularly for large r. 

• For a clean tube, 0 = Oo for water, 0 = 140o for mercury 

• For r > % in (6 mm), capillarity is negligible • Its effects are negligible in most engineering 
situations. 

• Important in problems involving capillary rise, e.g., soil water zone, water supply to plants 

• When small tubes are used for measuring properties, e.g., pressure, account must be made for 

ca P iUarit y Meniscus- 


c > 

t 




g | 

Water 


Mercury 


n 


4> 


(a) Wetting 
fluid 


(b) Nonwetting 
fluid 



v iy 







, * _ Meniscus- 
h> 0 

1 




V 1 




Water 



Mercury 

\ 

h < 0 
i 






[ 13 ! 

i ! 
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example ( 9 ) : If the surface tension of water-air interface is 0.069 N/m, what is the pressure 
inside the water droplet of diameter 0.009 mm? 

Solution: 


Given d= 0.009 mm; o= 0.069 N/m 

The water droplet has only one surface, hence, 


p = — = 

d 0.009X10 -3 


4X0.069 N rrn kN in 

= 30667 — = 30. 667 — or kPa 


m 


rrv- 


€xamole (W) : A clean tube of diameter 2.5 mm is immersed in a liquid with a coefficient of 
surface tension = 0.4 N/m. the angle of contact of the liquid with the clean glass can be assumed 
to be 135 0 . the density of the liquids 13600 kg/m 3 . what would be the level of the liquid in 
tube relative to free surface of the liquid inside the tube? 

Solution: 

Given d= 2.5 mm, o= 4 N/m, 0 = 135° ; p = 13600 kg/m 3 

Level of the liquid in the tube, h: 

2o" cos0 
h = — - 

pgR 

i 4 x 0.4 X cosl35 

h ~ (9.81 x 13600) x 2.5 x 10“ 3 
= — 3.3910 -3 m or — 3.39 mm 

Negative sign indicates that there is a capillary depression (fall) of 3.39 mm. 


14 
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mom Work : Derive an expression for the capillary height change h, as shown, for a fluid of 
surface tension a and contact angle 0 between two parallel plates W apart. Evaluate h for 
water at 20°C (g= 0.0728 N/m) if W - 0.5 mm. 


Solution: 

With b the width of the plates into the paper, the capillary forces on each wall together balance 
the weight of water held above the reservoir free surface: 

hbpgW = 2(abcos0 ) 

2(acos0) 
h PgW 


for water at 20°C (c= 0.0728 

u _ 2 x (0.0728 x cos(0)) 
h “ 9790 x 0.0005 


N/m, y — 
= 0.03m 


9790 N/m 3 ) and W = 0.5 mm. 
= 30mm 
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Chapter Wtot 

FLUID STATICS 


2.1. Introduction ; 

Fluid statics is that branch of mechanics of fluids that deals primarily with fluids at rest. 
Problems in fluid statics are much simpler than those associated with the motion of fluids, and 
exact analytical solutions are possible. Since individual elements of fluid do not move relative to 
one another, shear forces are not involved and all forces due to the pressure of the fluid are normal 
to the surfaces on which they act. Fluid statics may thus be extended to cover instances in which 
elements of the fluid do not move relative to one another even though the fluid as a whole may 
be moving. With no relative movement between the elements, the viscosity of the fluid is of no 
concern. 

In this chapter we shall first examine the variation of pressure throughout an expanse of fluid. 
We shall then study the forces caused by pressure on solid surfaces in contact with the fluid, and 
also the effects (such as buoyancy) of these forces in certain circumstances. 


2.2. Forces Exerted to an Element in a Fluid 


At a point, fluid at rest has the same pressure in all direction and to prove this, a small 
wedge-shaped free body element is taken at the point (x, y, z) in a fluid at rest. 
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16 i 



For unit width of element in z direction, and from the geometry of wedge we have the follows: 


Where 0 is an arbitrary angle, these results gives an important first principle of hydrostatics: 
At a point, fluid at rest has the same pressure in all direction. 

2.2.1. Pressure Variation: 

• For static fluid, pressure varies only with elevation (depth) change within fluid. 

To prove this real, we take a cubic fluid element as shown: 

While fluid at rest, applying the equations of equilibrium on the element. That’s yield: 

1. In vertical direction y 





Substitute of eq.3 in eqs. 1 and 2 and rearrange the terms yields: 



At a point the element limits to have an infinitesimal dimensions and then we can eliminate the 
1 

term (- dydx ■ y) from the above eq. because of it’s a higher order of differential values. Thus we 
have at final that: 




-> dP y ■ dxdz = —dxdydz ■ y 
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When the direction of Sy downward (nnj 17 In the negative direction of y ), which is 
called depth and denoted by h , last above" equation become: 

these results gives an important second principle of hydrostatics: 



• For static fluid, pressure varies only with elevation (depth) change within fluid by rate 
equal to specific weight y of that fluid. 

• In a fluid, pressure decreases linearly with increase in height 

• Pb - Pa = -y (zb - za) 

• This is the hydrostatic pressure change. With liquids we normally measure from the 
surface. 

• Open free surface pressure in liquids is atmospheric, Patmospheric. 

• For constant density fluids, and taking the free surface pressure as zero, p = y h. 

• Thus h = - 

y 

• Pressure related to the depth, h, of a fluid column. 

• Referred to as the pressure head, h. 
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2. In horizontal direction x : 



E f r = Px ' dydz - (T> + dPx) ■ dydz = 0 


P r ■ dydz — P r ■ dydz — dP r • dydz = 0 .... (2. 9) 


-> dPy = 0 ... . (2. 10) 


This equation means there is no change in horizontal pressure with horizontal direction. 


Example i : A freshwater lake near Bristol, New Hampshire, has a maximum depth of 60m, 
and the mean atmospheric pressure is 91 kPa. Estimate the absolute pressure in kPa at this 
maximum depth. 

Solution: 

From Table 2.1, take y = 9790 N/m3. With P a = 91 kPa and z = - 60 m, the pressure at this depth 
will be : 

P = 91 kN/m 3 - (9790 N/m 3 )(-60 m ) 

v j\ j 1000 N 

P = 91 kPa + 587 kN/m 2 = 678 kP a 

By omitting P a we could state the result as P = 587 kP a (gage). 


@xatnvle 2: The liquid oxygen(LOX) tank of a Saturn moon rocket is partially filled to a depth 
of 30 ft with LOX at - 320 °F . The Pressure in the vapor above the liquid surface is maintained 
at 14.7 psia. Calculate the pressure at the inlet valve at the bottom of the tank. 

Solution: 

y = p * g = 2.33 *32.2 = 75 lb /ft 3 

Pi = 75 * 30 + 14.7 * 144 = 2341 psfa = 16.3 psia 

Calculate pressure and specific weight at air in the U.S. standard atmosphere at altitude 35,000 
ft. 
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Solution: 

At sea level : 

Pi = 14.7 psia , Ti = 519° R , yi = 0.0765 lb/ft 3 

P _ 14 . 7*144 1 _ 6450 

yi.235 — ( 0 . 0765) 1 - 235 ' 7 ~~ P 0 - 81 

| 19 | 

35,000 = 6450 J 14 - 7 * 144 p- o.ei dp ? p 2 = 504 ' ponr -i 3.5 psia 

T 2 = 519 - 35,000 * 0.00356 = 394° R 
y 2 = 504 * 32.2 / 1715 * 394 = 0.024 lb/ft 3 


2.2.2. Absolute and Gauge Pressure : 

In practice, pressure is always measured by the determination of a pressure difference. If the 
difference is that between the pressure of the fluid in question and that of a vacuum then the result 
is known as the absolute pressure of the fluid. More usually, however, the difference determined 
is that between the pressure of the fluid concerned and the pressure of the surrounding 
atmosphere. 

This is the difference normally recorded by pressure gauges and so is known as gauge pressure. 
•If P < Patm , we call it a vacuum (or negative or suction) pressure , 
its gage value = how much below atmospheric 

•Absolute pressure values are all positive 

•While gage pressures may be either : 

- Positive: if above atmospheric, or 

- Negative (vacuum, suction): if below atmospheric 
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•Relationship between absolute, gage and atmospheric pressure reading: P a bs = Patm + Pgage 



Pascal (Pa) 

bar (bar) 

technical 

atmosphere 

(at) 

atmosphere 

(atm) 

torr (Torr) 

pound-force /in 2 

(psi) 

1 Pa 

= 1 N/m 2 

10“ 5 

1.0197xl0“ 5 

9.8692xl0“ 6 

7.5006xl0“ 3 

145.04X10 -6 

1 bar 

100,000 

= 10 6 dyn/cm 2 

1.0197 

0.98692 

750.06 

14.5037744 

1 at 

98,066.5 

0.980665 

= 1 kgf/cm 2 

0.96784 

735.56 

14.223 

1 atm 

101,325 

1.01325 

1.0332 

= 1 atm 

760 

14.696 

1 torr 

133.322 

1.3332xl0 -3 

1.3595xl0 -3 

1.3158xlO -3 

= 1 Torr; 
~ 1 mmHg 

19.337xl0 -3 

1 psi 

6.894xl0 3 

68.948xl0 -3 

70.307xl0 -3 

68.046xl0" 3 

51.715 

= 1 M/in 2 


example 3: 1 Pa = 1 N/m2 = 10-5 bar = 10.197x10-6 at = 9.8692x10-6 atm > etc. 


% 

$ 


Atmospheric pressure 

Vacuum = negative 
gage pressure 

Absolute 

pressure 

Absolute zero | 


Gage 

pressure 


Absolute 

pressure 


Atmospheric 

pressure 


Pressure Units 
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A pressure of 1 atm can also be stated as: 

= 1.013 25 bar 

= 1013.25 hectopascal (hPa) 

= 1013.25 millibars (mbar, also mb) 

= 760 torr [B] 

~ 760.001 mm-Hg, 0 °C~ 1.033 227 452 799 886 kgf/cm 2 
~ 1.033 227 452 799 886 technical atmosphere 
~ 1033.227 452 799 886 cm-H 2 0, 4 °C 

I 21 ! 

I 1 

€xamule 4: For two points A, B: PA - 200 kPa gauge , PB - 51 kPa abs . , where P a tm = 101 
kPa , Re-read the value of pressures in A and B in all possible forms. 

? A - 200 kPa gauge 

Pa=200 kPa gauge 


P atm. - 0 kPa 
P b =51 kPa 
Pabs = 0 kPa 


P a =301 kPa abs. 


P B = -50 kPa gauge 
P b =50 kPa vacuum 
Pb= 50 kPa suction 


P b =51 kPa abs. 


Example 5 : In Fig. the tank contains water and immiscible oil at 20°C. What is h in cm if the 
density of the oil is 898 kg/m3? 


"D%. "padAil /iid dfl- 


"plaid "7KecA<zmc4 


Secmd Stacie 


Solution: For water take the density - 

998 kg/m . Apply the hydrostatic relation 
from the oil surface to the water surface, 
skipping the 8-cm part: 

Patm + (898)(g)(h + 0. 12) 

-(998)(g)(0.06 + 0.l2) = p atm . 

Solve for h ~ 0.08 m = 8.0 cm Ahs. 


h 

12 cm 


6 cm 


Oil 



8 cm 


Water 
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ffiorne Worbg : 

1. A closed circular tank filled with water and connected by a U-piezometric tube as shown 
in figure. At the beginning the pressure above the water table in the tank is atmospheric, 
then the gauge that connected with tank read an increasing in pressure that caused falling 
in the water level in the tank by 3 cm. a) calculate the deference in height that accrued 
between water levels inside the tank and in the external tube leg. b) Determine the final 
pressure that was reading by the gauge. 


2. One means of determining the surface level of oil in a tank is by discharging a small amount 
of air through a small tube, the end of which is submerged in the tank, and reading the 
pressure on the gauge that is tapped into the tube. Then the level of the oil surface in the 
tank can be calculated. If the pressure on the gauge is 30 kPa and the specific gravity of oil 
is S = 0.85, what is the depth d of oil in the tank? 


23 
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2. 3. Pressure measurement devices : 

2.3.1. Absolute pressure measurement (Barometers) : 

Barometers : The instrument used to measure atmospheric pressure is called barometer 

1. Mercury Barometer : which is illustrated in figure below, which consist of a one meter 
length tube filled with mercury and inverted into a pan that’s filled partially with 
mercury. The height difference of mercury in inverted tube respect to outside them reads 
the atmospheric pressure value. 

- Values of standard sea-level atmospheric pres sure=7 01. 325 kPa abs =1013.25 mbar abs 
= 760 mm Hg 
=10.34 m H 2 0 



(a) 


(b) 


A barometer measures local absolute atmospheric pressure: ( a ) the height of a mercury column is proportional 
to /Mm; (b) a modern portable barometer, with digital readout, uses the resonating silicon element 
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2. Aneroid barometer: uses elastic diaphrj 24 j> measure atmospheric pressure 

! ■ 





(Example 6: (Piezometer) 

In figure pressure gage A reads 1.5 kPa. The fluids are at 20oC. Determine the elevations z, 
in meters, of the liquid levels in the open piezometer tubes B and C. 


V 

I ... 

L 



Solution : (B) Let piezometer tube B be an arbitrary distance H above the gasoline- glycerin 
interface. The specific weight are y a ir = 12 N/m 3 , y ga soiine = 6670 N/m 3 , and y g i yC erin= 12360 
N/m 3 . then apply the hydrostatic formula from point A to point B. 

1500 N/m 2 + (1 2.0 N/m 3 )(2.0 m) + 6670( 1 .5 - H) - 6670(Z B - H - 1 .0) = p B =0 (gage) 

Solve for Z ft = 2.73 m (23 cm above the gasoline-air interface) Ans. (b) 

Solution (C): Let piezometer tube C be an arbitrary distance Y above the bottom. Then 

1500 + 12.0(2.0) + 6670( l .5) + 12360(1.0 - Y) - 12360(Zc - Y) = p c = 0 (gage) 

Solve for Z c = 1.93 m (93 cm above the gasoline-glycerin interface) A ns. (c) 


Air 


Gasoline 


Glycerin 


' 
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2.3.2. Gage pressure measurement: 

2.3.2. 1. Manometry 

1. Piezometer 

For measuring pressure inside a vessel or pipe in which liquid is there, a tube open at the top 
to atmosphere may be attached, tapped, to the walls of the container (or pipe or vessel) containing 
liquid at a pressure (higher than atmospheric) to be measured, so liquid can rise in the tube. By 
determining the height to which liquid rises and using the relation PI = pgh, gauge pressure of 
the liquid can be determined. Such a device is known as piezometer. To avoid capillary effects, a 
piezometer's tube should be about 12mm or greater. 

Open 



2. Manometers 

For high pressures measurement using piezmeters become impractical and it is useless for 
pressure measurement in gases. The manometers in its various forms is an extremely useful type 
of pressure measuring instrument for these cases. 


Open 



7 

h 2 

1 


( 3 ) 
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Professor John Foss (Michigan State University) Procedure for manometers pressure 
calculation: 


. Manometers limitations: manometers suffers from a number of limitations. 

1. While it can be adapted to measure very small pressure differences, it cannot be used 
conveniently for large pressure differences - although it is possible to connect a number of 
manometers in series and to use mercury as the manometric fluid to improve the range, 
(limitation) 

2. A manometer does not have to be calibrated against any standard; the pressure difference 
can be calculated from second and third principles in hydrostatics. ( Advantage) 

3. Some liquids are unsuitable for use because they do not form well-defined interface. 
Surface tension can also cause errors due to capillary rise; this can be avoided if the 
diameters of the tubes are sufficiently large - preferably not less than 12 mm diameter, 
(limitation) 

4. A major disadvantage of the manometer is its slow response, which makes it unsuitable for 
measuring fluctuating pressures 

5. It is essential that the pipes connecting the manometer to the pipe or vessel containing the 
liquid under pressure should be filled with this liquid and there should be no air bubbles in 
the liquid. (limitation). 
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Curved tube of elliptical cross-section changes curvature with changes in pressure. 
Moving end of tube rotates a hand on a dial through a linkage system. Pressure indicated 
by gage graduated in kPa or kg/cm 2 (=98.0665 kPa) or psi or other pressure units. 



example 7 : (U -manometers) 

In Figure fluid 1 is oil (Sg=0.87) and fluid 2 is glycerin at 20oC (y=12360 N/m3). If Fa-98 
kPa, determine the absolute pressure at point A 



I 
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exzmvle 8: (Differential-Manometers)Pressure gage B is to measure at point A in a water flow. 
If the pressure at B is 87 Kpa, estimate the pressure at A, in Kpa. Assume all fluids are at 20° 


C. 



Solution 


First list the specific weights from Table 2. 1 or Table A, 3: 

Twaler = 9790 N/lll' Tmcrcury “ 133,100 N/llV 7 0 j| = 8720 NW 
Now proceed from A to 5, calculating the pressure change in each fluid and adding: 

Pa ~ 7w(te)w ~ Y ~ = Pb 

or p A - (9790 N/m 3 )(- 0.05 m) - (133,100 N/m :, )(0.07 m) - (8720 N/m 3 )(0.06 m) 

- p A + 489-5 Pa - 9317 Pa - 523.2 Pa - p B = 87,000 Pa 

where we replace N/m“ by its short name. Pa. The value — 0.07 m is the net elevation 
change in the mercury (1 1 cm — 4 cm). Solving for the pressure at point A, we obtain 

p A = 96,35 1 Pa = 96.4 k Pa Ans. 
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In the figure below, all fluids are at 20°C. Determine the pressure difference (P a ) between 
points A and B. 



Solution: Take the specific weights to be 

Benzene: 8640 N/m^ Mercury: 133100 N/m 3 

Kerosene: 7885 N/nV* Water: 9790 N/m^ 

and % ir will be small, probably around 12 N/nV . Work your way around from A to B: 

p A + (8640)(0.20 m) -(1331 00x0.08) - (7885)(0.32) + (9790)(0.26) - (1 2)(0.09) 
— p B , or, after cleaning up, p A — p B - 8900 Pa Arts. 


€xamole 10: ( Inverted-Manometers ) 

For the inverted manometer of the figure below, all fluids are at 20°C. If Pb - Pa = 97 kPa, 
what must the height H be in cm? 



i 
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Solution: Gamma — 9790 N/m' for water 

and 1 33 L 00 N/m for mercury and 
(0.827X9790) - 8096 N/m^ for Mcriam red 
oil. Work your way around from point A to 
point B : 

p A -(9790 N/m 3 )(H meters) -8096(0. 18) 

+ 133 100(0. 18 + 11+ 0.35) = p B = p A +97000. 

Solve tor tf «= 0.226 m = 22.6 cm A ns. 


€xanwle 11: (inclined- Man o meters (micro-manometer) 

The ratio of the cistern diameter to tube diameter is 10, when air in the cistern is atmospheric 
pressure, the free surface in the tube is at position 1. When the cistern is pressurized, the liquid 
in the tube moves 60 cm up the tube from position 1 to position 2. What the cistern pressure 
that cause this deflection? The specific gravity of the liquid is 1.2. 



Solution: 


P 1 = P 2 + 1.2 x 9810 x 0.6 x sin (30) = 11772 x 0.3 = 3531.6 Pa 
P 'cistern = P\ + h X 1.2 X 9810 

P cistern 10 ^ tube 

cistern = tube 

D 2 d 2 

X 7T X h = X 7T X 0.6 

4 4 

100 x d 2 d 2 

X 7T X h = X 7T x 0.6 

4 4 

h = .006 = 0.6 cm 

Pcistern = 3531.6 + 0.006 X 1.2 x 9810 = 3602.2 Pa 
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€xamvle 12: (U -Manometer) 

This vertical pipe line with attached gage and manometer contains oil and mercury as 
shown. The manometer is open to the atmosphere. There is no flow in the pipe. What will 
be the gage reading P x ? 



Solution : 

Pi = P x + (0.9 * 62.4) * 10 
Pi = (13.57 *62.4) *(15/12) 
Because 
Pl = Pr 

P x = 497 psf = 3.45 psi 
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Chapter Qfyxtt 

HYDROSTATIC forces 


3.1. Introduction ; 

The pressure of a fluid causes a thrust to be exerted on every part of any surface with which 
the fluid is in contact. The individual forces distributed over the area have a resultant, and 
determination of the magnitude, direction and position of this resultant force is frequently 
important. For a plane horizontal surface at which the fluid is in equilibrium the matter is simple: 
the pressure does not vary over the plane and the total force is given by the product of the pressure 
and the area. Its direction is perpendicular to the plane - downwards on the upper face, upwards 
on the lower face - and its position is at the centroid of the plane. But if the surface is not 
horizontal the pressure varies from one point of the surface to another and the calculation of the 
total thrust is a little less simple. 

3. 2. Hydrostatic Forces on Fianc Surfaces : 

3.2.2. 1. On Horizontal Surface: 

Free surface 


P = Palm \ v 


i 

i 

Specific weight = y 

/ 


Fr 



P 

\ 

j \ \ \ \ ' 

' \ \ 1 \ \ I \ 


~P Palm f\ 


Resultant force is : 


Pp = PA = yh A ....( 3 . 1 ) 


A: the bottom area of container 
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3.2.2.2. On an Inclined Surface: 



Hydrostatic force and center of pressure on an arbitrary plane 
surface of area A inclined at an angle 0 below the free surface. 

Consider a plane shown above : 

- At surface: p - p a m 

- Angle 6 between free surface & the inclined plane 
• Along the vertical depth h 

^ Pressure linearly changes ^ Hydrostatic force changes 
Differential Force acting on the differential area dA of plane, (Perpendicular to plane): 


dF =(pressure).(Area)=(yh).(dA) .... (3.2) 


Then, Magnitude of total resultant force Fp 
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Fp=Ja2_ ( sin e _ ) dA - ■ ( 3 - 3 ) 


\F P = y sin 9 f. y dA ... . (3.4) 


L yd . a 

y c = A - 


f, x dA 

V — A 

ZE a " 

••(3.6) 


where h = y sin 6 

(1st mi" ! of the area about x-axis ) 


y c , x c : y and x coordinates of the center of area respectively (Centroid) 


MY C = f ydm ....(3.8) 


Yc.Xc: y and x coordinates of the center of mass respectively (first moment) 


1 = L r2dm ■■■ ■ (3.9) 


( second moment of mass) 


i x = L y 2dA ••••(3.10) 


I y = f x 2 dA ....(3.11) 


l x , Iy : second moment ( moment of inertia ) about x and y coordinates respectively, then: 


Pp = y A yc sin 9 = ( y he ) A .... (3.12) 


where y h c : Pressure at the centroid 


y h c = (Pressure at the centroid) x Area 
The location of point of action of Fr 

- Not passing though Centroid & related with the balance of torques due to of Fp 
i) Position of Fp on y-axis 

yp : y coordinate of the point of action of Fp 
Taking moment about x axis'. 
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The moment of resultant force = The moment of its components 

F P y P = I ydF 

\ 35 ; 

{/Ay c s jft0)y R = \jffn0y 2 dA = y^nhyj-'dA 



f 2 

where I x = y dA: 2nd moment of area (Moment of inertia, +ve always ) 

2 

or, by using the parallel-axis theorem, I x = I xc + Ay c 

(Always below the centroid ) 

ii) Position of FP on x-axis 

xp : x coordinate of the point of action of Fp 

The moment of resultant force = The moment of its components 

F p x p = \xdF 

1 v> i/Ay c ^fn 6)x p = yytn OxydA - y /m 0j A xydA 



where J A xydA = I xy : Area product of inertia (+ve or -ve) 
or, by using the perpendicular-axis theorem, I xy ~ ^ xyc + Ax c y c 


Xp = x c ....(3.16) 

yrA 


c.f. For a symmetric submerged area, xp = x c (Ixyc = 0) 
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Moments of Inertia of Several Common Cross Sections 
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@xattwle l :_ The gate in Fig. below is 5 ft wide, is hinged at point B, and rests against a smooth 
wall at point A. Compute (a) the force on the gate due to seawater pressure, (b) the horizontal 
force P exerted by the wall at point A, and (c) the reactions at the hinge B. 


Wall 


Solution : 



Part a : 

By geometry the gate is 10 ft long from A to B, and its centroid is halfway between, or at 
elevation 3 ft above point B. The depth hcG is thus 15- 3 = 12 ft. The gate area is 5(10) = 50 ft 2 . 
Neglect p a as acting on both sides of the gate. From Eq. (3.4) the hydrostatic force on the gate is 


F = pcoA = yhccA = (64 lbf/ft 3 )(12 ft)(50 ft 2 ) = 38,400 lbf Ans. (a) 

Part b : 

First we must find the center of pressure of F. A free-body diagram of the gate is shown in 
figure below. The gate is a rectangle, hence 
I xy = 0 and Ixx = bL 3 /12 = 5*10 3 /12 = 417 ft 4 

The distance / from the CG to the CP is given by Eq. (3.10) since p a is neglected. 

/ = - yep = + 


h cc A 
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The distance from point B to force F is thus 10-/-5 = 4.583 ft. Summing moments 
counterclockwise about B gives : 

PL sin 0 - F(5-/) = P*6 - 38400 * 4.583 = 0 
P = 29300 ft. 

Ans. (b) 

Part c : 

With F and P known, the reactions Bx and Bz are found by summing forces on the gate 
S F x = 0 = B x + F sin 0 - P = B x + 38400 * 0.6 - 29300 
B x = 6300 lbf 

I F z = 0 = B z - F cos 0 = B z - 38400 * 0.8 

B z = 30700 lbf Ans. (c) 


€xamole 2: A tank of oil has a right-triangular panel near the bottom, as in figure below Omitting 
p a , find the (a) hydrostatic force and (b) CP on the panel. 


Solution : 
Part a : 



The triangle has properties given above. The centroid is one-third up (4 m) and one-third 
over (2 m) from the lower left corner, as shown. The area is 
A= 0.5*6*12 = 36 m 2 
The moment of inertia are : 

! 39 ! 


*D%. 'padAil /f&d /ft- /ftf&ao ^facccut 


*ptuid ‘TftecAanico 


Second Stooge 


Ixx = bV/36 = 6*12 3 /36 = 288 m 4 


W 


b(b-2s)L 2 


6 [ 6 — 2 * 6 ]* 12 2 


= -72 m 4 


72 72 

The depth to the centroid is he g = 5 + 4 = 9 m; thus the hydrostatic force from Eq. (3.10) is 
F = p g hco A = 800* 9.807 * 9 * 36 = 2.54* 10 6 (kg.m)/s 2 = 2.54* 10 6 N = 2.54 MN Ans. (a) 


Part (b): 


yep = - 


h CG A 


288 sin 30° 
9*36 


= —0.444 m 


x C p = 


Ixy sin 0 


(-72) sin 30° 


= +0 111 m 


Ans. (b) 


hcc;A 9*36 

The resultant force F = 2.54 MN acts through this point, which is down and to the right of the 
centroid, as shown in figure above. 


Example 3: Calculate magnitude, direction and location of the total force exerted by the water on 
one side of this composite area which lies in a vertical plane. 
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Solution: 

By inspection the direction of the force is normal to the area. 

Magnitude : 

Force on triangle = 62.4*3*12 = 2245 lb 
Force on quadrant = 62.4*3. 55*(367i/4) = 6275 lb 
Total force of composite area = 2245 + 6275 = 8520 lb 
Vertical location of resultant force : 

— for triangle = — - — — = 0.667 ft 

l c A 36 * 3*12 

/ 711 

— for quadrant = : — = 0.71 ft , taking moment about line AC, 

3 . 55 * 28.3 

2245*2.667 + 6275*3.26 = 8520 (Z p - 1) 

Z P = 4.09 ft 

lateral location of resultant force : since the center of pressure of the triangle is on the median line 
, x p i , is given (from similar triangles) by 

( 1 . 333 - x pl _ 0.667 

2 ~ 6 
Xpl = 1 . 1 1 1 ft 

Dividing the quadrant into horizontal strips of differential height, dy , the moment about BD of 
the force on any one of them is : 
dM = (x dy ) y h (| ) 

in which h = y + 1 and x 2 * * * + y 2 = 36, substituting and integrating gives the moment about BD of 
the total force on the quadrant, 

M = ^ / 0 6 ( 36 - y 2 )(y + l)dy = 14600 ft. lb 
and thus 

x p2 = 14600/6275 = 2.33 ft 
finally taking moment about line BD, 

2.33 * 6275-2245*1.111 = 8520 * x p 
x p = 1 .425 ft 

thus the center of pressure of the composite figure is 1 .425 ft to the right of BD and 4. 1 ft below 
the water surface. 


I 


41 


*D%. 'padAil /f&d /ft- /ftt&ao ^facccut 


*ptuid ‘TftecAanico 


Second Stacie 


example 4: An inclined, circular gate with water on one side as shown in figure. Determine the 
total resultant force acting on the gate and the location of the centre of pressure. 

Solution: 


F = y h cg A = (9.79)[1.5+0.5 (tt( 1) 2 /4] = 14.86 kN 


Z C p — Zc2 + 


'Cg 


l cg 

%cgA 


=[- 


60 


+ 0.5 * 1 + 


7T 


(l ) 4 

64 


[1.5/ sin 60+0.5*(1 )][tt*^H 

4 


2.26 m 


60 * 



3. 3. Hydrostatic Forces of Pressure Prism : 


Especially applied for a rectangular surfaces (areas) , a Simple method for finding the force and 
the point of action. 


Consider the situation shown : 



JL 
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Information from the diagram 

- Vertical wall of width b and height h 

- Contained liquid with specific weight y 

- Pressure. Pf Q p ^ ^ Pbottom — Y h 


From the last section, 


F p = ( yh c ) ■ (A) = p av (at the centroid) x area = 



A 


Let’s define a pressure-area space. (See the right figure above] 

1. Horizontal axis: Magnitude of the pressure 

2. Vertical axis: Height of the area 

3. Axis toward the plane: Width of the area 

: Resultant volume ( Pressure prism) 

• How to find the resultant force Fr from the pressure prism 


Fp = lyhA = iyh(bh)....(3.17) 


(Volume of the pressure prism) 


3.3.1. Finding the point of action ofF/t (the point of action) : 


From the last section , equation 3.14 : 


Ac i Ac 


y P = — + y r = 

VrA 0.5ft(bfa) 


+ 0.5ft 


- In case of rectangular plate : 
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Then eq. 3.14 becomes : 
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' = $S + i fc= ; fc+ 5 ft= i fc ---( 3 - 19 > 


(from the top) 


From the pressure prism : 

Yp = vertical center ( Centroid of the pressure prism ) 

2 , 1 , 

= —n (from the top) = — n (above the base) 

3 3 

Xp = Horizontal center 



Note: Special case of a plane surface not extending up to the fluid surface 

- Completely submerged plane . Consider the situation shown : 

1. Resultant force Fp 

Fp = Volume of the shadow region 
Fp = Volume of hexahedron + Volume of prism 


Fp = Fi(abde) + F 2(bcd) .... ( 3 . 20 ) 


7 p =(yhi)A + 0.5[y (h 2 - hi)] A .... ( 3 . 21 ) 



y(/j 2 - hj) 


2. The location of Fr (y a) : Consider the moments again 

Moment by F p acting at y A = Moment by F [ at + Moment by at y 2 

h 2 h 

F p y p = F ] y, + F 2 y 2 where y\=— (for rectangle) y 2 = — (for triangle) (From the top) 

2 3 

• The effect of the atmospheric pressure p atm 

: Increasing Volume of hexahedron (F\), NOT the prism (Fi) 
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exzmvle 5 : The dam of figure has a strut AB every 6 m. Determine the compressive force in the 
strut, neglecting the weight of the dam. 




3 tn — ► 


Solution: If the struts are doing their job. the moment about the hinge should be 0 ! Every six 
meters there is a strut, so the width eaeh strut supports is 6 m. 


Use any method to determine the forees on the dam due to the water. Using the pressure prism 
method: 

( 1 ) the foree on the upper 2 m portion of the dam is 

F 1 — ^ (27) (2 m) (6 m wide) — 12 7 
(2) the foree on the lower portion of the dam is split into 2 pieees 

t'l — (27) (5 m) (6 m wide) — 60 7 

F-i — 7j (47) (5 m) (6 m wide) = 60 7 


and 
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[Tic compressive force on the strut AB is in a direction along the strut and so to determine the 
moment about the hinge, the force is split into x- and y-eomponents. determined by the geometry 
shown in Fig. 2.6 h 

p V§0 „ .,4 

tf AB x ~ 


6 


AB 


*AB y = - F A B 


Now writing the moment equation about the hinge, (moment = force x lever ann) 


+ 0 EMu - 0 : 


-Fi 


4 + |( 2 ) 


F 2 (2.5) -ft 


5 ® 


+ - i |^^AB(4) + ^AB(3)-C 


or, alter rearranging 


4.9814 F A b = 549, 186 N-m + 1,470,900 N-m + 980,600 N-m 
^ab = 602,368 N - 602.4 kN 


(Example 6 : A structure is so arranged 45 Jg a channel that it will spill the water 

out if a certain height y (Fig. 2.18a) is lied. The gate is made of steel plate 

weighing 2500 N/m 2 . Determine the height of y. 



SOLUTION Using pressure-prism concepts, for unit width normal to the page the 
force on the horizontal leaf (Fig. 2. 1 86) is given by the volume of a pressure prism 
of base 1.2 m 2 and constant altitude yy N/m 2 , which yields F y — \.‘ly y N acting 
through the center of the base. The press tiro prism for the vertical face (Fig. 2. 18c) is 
a wedge of base y m 2 and altitude varying from 0 to y y N/m 2 . The average altitude 
is 7^/2, so F x = yy 2 /‘2 N. The centroid of the wedge prism is y / 3 from the hinge. 
The weight of the gate floor exerts a force of 3000 N at its center. Figure 2. 1 %d sliows 
all tiie forces and moment arms. For equilibrium, that is, the value of y for tipping, 
moments about the hinge must be zero. 
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or 


M = (3000 N)(0.6 m) + (1 :2yy N)(0.6 m) 



M - y s -4:S2y- 1.1014-0 

By using try and error technique (or other techniques like Newton Raphson method) we find that 
y=2.196m 


46 
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3. 4. Hydrostatic Forces on Carved Surfaces : 
Case no.l: 





For unit width of surface 



dF P y = P ds cos 6 
P = yy and dx = ds cos 0 
dFpy = yy dx 


:.i F py = f yy dx = y / y dx = y V ... . (3.23) 


Where V the volume of liquid above the surface to the zero pressure surface 


I 


47 
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Case no.2: 

By the same way we find the vertical component to pressure force if the liquid exist lower the 
surface by taking the sign of V as -ve to represent the upward direction of this force 




dFpy = P ds cos 0 = - P ds cos 0 
P = y y & dx = ds cos 6 

dFpy = - yy dx 


f yy dx = —yfydx = —yV.... (3.24) 


dFpx = P ds sin 0 
P = yy & dy = ds sin 6 
dF px = yydy 

• w f px = Sy* Yydy = y f** y dy 


y y 4 ..-.(3.25) 


r * 


yy z 

Where ( is the volume of pressure prism on the surface projection on vertical plan 



The Magnitude of Fp is: 



in direction of 


tan 9 = — 

Fpx 


& the line of action can be find from the concept of: 

Moment of resultant force = The summation of the moments of its components, i.e. 
MF r =Y d (MF r ,,MF„) 
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example 7: A curved surface AB is a circular arc in its section with radius of 2m and width of 
lm into the paper. The distance EB is 4m. The fluid above surface AB is water, and 
atmospheric pressure applied on free surface of water and on the bottom side of surface AB. 
Find the magnitude and line of action of the hydrostatic force acting on surface AB. 


Solution: 



1 . Equilibrium in the horizontal direction 

F x = F h = P A = (5 m )* (9810 N/m 3 )*(2*1 m 2 ) = 98.1 kN 


2. Equilibrium in the horizontal direction 

■ Vertical force on side CB 

F,. = pyA = 9.81 kN l m 3 x 4 m x 2 m *1 m. = 78 .5 kN 

■ Weight of die water in volume ABC 

— (9 81 kl'I/m 5 ) x (0 25 x t x 4 m 3 ) ( l m) = 30 8 kN 

■ Summing forces 

F r = W + Py= 109 3 kN 

3. Line of action (horizontal force) 

Vcp = 5.067 m 

4. The line of action (x C p) for the vertical force is found by summing moments about point C: 

HcpF y — F y x 1 m 4 W xxjp 

The horizontal distance from point C to die centroid of the area ABC is found using Fig. A. 1 : x 
w = 4r/3D= 0.849 m. Thus, 

„ _ 73 5 kN* 1 m + 30.8 kN x 0 S49 m _ ^ 

S v 109 1 kN V - 9j,a 

49 ! 


i 
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5. The resultant force that acts on the curved surface is shown in the following figure. 



example 8: A cylindrical barrier in Fig. holds water as shown. The contact between the cylinder 
and wall is smooth. Consider a 1-m length of cylinder; determine (a) its weight, and (b) the 
force exerted against the wall. 



SOLUTION (<x) For equilibrium the weight of the cylinder must equal the vertical 
component of force exerted on it by the water. (The imaginary free surface for CD is 
at elevation -4). The vertical force on BCD is 


Fvucu = + 2r2 ) T = ( 27r + 8 ) T 


The vertical force on AB is 


F, 


VAIS 



.'ll:? i| 

50 ! 


7 = — (4 — 7r) 7 
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Hfciice, tiie weight per meter of length is 

F vacn + F Vau (3tt 4- 4)7 - 0.132 MN 

(ft) The force exerted against the wail is the horizontal force on ABC minus the horizontal 
force on CD. The horizontal components of force on BC and CD cancel; the projection 
of BCD on a vertical plane is zero. Hence, 

F h = F HAa = 2 7 = 19.6 kN 

since the projected area is 2 m 2 and tiie pressure at the centroid of the projected area is 
9806 Pa. 


51 ! 
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Chapter Jfour 

ACCELERATED FLUIDS WITH 
LINEAR AND ANGULAR ACCELERATION 

4./. Introduction : 

If a body of fluid is moved at a constant velocity, then it obeys the equations derived earlier for 
static equilibrium. 



If a body of fluid is accelerated such that, after some time, it has adjusted so that there are no 
shearing forces, there is no motion between fluid particles, and it moves as a solid block, then the 
pressure distribution within the fluid can be described by equations similar to those applying to 
static fluids. 

Consider the forces acting on a small horizontal element, area 8A and length §x, with a uniform 
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acceleration ax in the x direction: 


F x = p8A- 



= - -T 8x8A....(4.2) 

dx v y 


But Newton's 2nd law gives: 


Px = m a x 



Now looking at the forces acting on a small vertical element, area SA and length Sy, with a uniform 
acceleration az in the z direction: 


F y = [ p 8A — (p + dy ) & 4 J — p 8y 8A g = p 8y 8Aa y ... . (4.5) 


And hence : 


j^= ~ P (.9 + a y) ■■■ ■ (4-7) 


Here, pressure p, is a function of x and y, so, by definition: 


\dp = ^ 8x 4- ^ 8y 
dx dv 


(4.8) 


and identifying isobars as lines of constant pressure with dp = 0, we have an equation for the 
slope of an isobar as: 


dp = dp dy f49 . 

dx dyjdx _ ' _ 


Then : 
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dp 



dy dx 

P a x 

..( 4 . 10 ) 

dx ^2. 

dy 

p(a+a y ) 


It is clear that if a x -J 0, then the isobars will not be horizontal in this case. 


L.2. Uniform Linear Acceleration 


i 53 1 

If a container of fluid is accelerated uniformly ; ^rizontally with a y = 0, then the slope of the 

isobars within the fluid is given by: 



Notes : 

• hydrostatic pressure variation in a gas is negligible in comparison to that in a liquid, 

• the free surface of a liquid is normally taken as a line of constant pressure - or isobar - and 


the equation above gives the surface slope of an accelerated fluid. 
• The resultant acceleration perpendicular on the isobar surfaces is: 



€xamole (1} A drag racer rests her coffee mug on a horizontal tray while she accelerate at 7 
m/s 2. The mug is 10 cm deep and 6 cm in diameter and contains coffee 7 cm deep at rest, (a) 
Assuming rigid body acceleration of the coffee, determine whether it will spill out of the 
mug.(b) calculate the gage pressure in the corner at point A if the density of coffee is 1010 
kg/m3. 
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Solution : 

The free surface tilts at the angle 0 gii 54 j eq. 4. 1 1 regardless of the shape of the mug. 
With a y = 0 and standard gravity : 1 

9 = tan~ 1 —= tan -1 — =35.5° 

g 9.81 

If the mug is symmetric about its central axis, the volume of coffee is conserved if the tilted 
surface exactly at the centerline, as shown in figure above : 

Thus the deflection at the left side of the mug is : 

y = 3 * tan 0 = 2.14 cm Arts, (a) 

this is less than the 3 cm clearance available, so the coffee will not spill unless it was sloshed 
during the start- up of acceleration. 

When at the rest, the gage pressure at point A is given by : 

Pa = P g (y surface - y A ) = 1010 * 9.81* 0.07 = 694 N/m 2 = 694 Pa. 

During acceleration, eq. 4.12 applies, with : 

G = [(7) 2 + (9.8 1) 2 ] 05 = 12.05 m/s 2 

The distance As down the normal from the tilted surface to point A is : 

As = ( 7 + 2.14 ) cos 0 = 7.44 cm. 

Thus the pressure at point A becomes: 

Pa = pG As = 1010 * 12.05 * 0.0744 = 906 Pa Ans. (b) 

Which is an increase of 31% over the pressure when at rest. 


€xattwle (2) : An open tank of water is accelerated vertically upward at 15 ft/sec 2 . Calculate the 
pressure at a depth of 5 ft. 


Solution : 


*D%. 'padAil /f&d /ft- /fO&ao ^facccui 


*ptuid ‘TftecAanico 


Second Stooge 


dp 

dz 


-62.4 ( 15 3 + 2 3 2 2 - 2 ) = -91.5 lb/ft 3 


P = C d P = “ f 0 5 91 - 5 dz = 457 P sf 


example (3) : This open tank moves up the\ 55 \ with constant acceleration .Calculate 
acceleration required for the water surface To iriove to the position indicated. Calculate 
pressure in the corner of the tank at A before and after acceleration. 



the 

the 


From geometry, the slope of the water surface during acceleration is -0.229. from the slope of the 
plane, a x = 4 a z : 

—0.229 = — — — 

(a z +32.2) 

And from the foregoing : 

a z = 1.96 ft/sec 2 , a x = 7.84 ft/sec 2 , a = 8.08 ft/sec 2 

from geometry, the depth of water vertically above comer A before acceleration is 2.91 ft ; hence 
the pressure there is 2.91 * 62.4 = 181.5 psf. After acceleration this depth is 2.75 ft : 

t = -( 32 - 2 + !- 96 ) (ii) = - 66 -° 

Therefore : 


P = 2.75 *66 = 181.5 psf. 

Note: 


The fact that the pressure at point A are the same before and after acceleration is no coincidence; 
general proof may be offered that pa does not change whatever the acceleration. This means that 
the force exerted by the end of the tank on the water is constant for all accelerations. However, 
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this is no violation of Newton’s second law , since the mass of liquid diminishes with increased 
acceleration so that the product of mass and acceleration remains constant and equal to the applied 
force. 


4.3. Uniform Angular Acceleration^^ 56 j w about Vertical Axis): 


If liquid is placed in a container and rotated about a vertical axis at a constant angular velocity, 
then after some time it will move as a solid body with no shearing of the fluid. Under such 
conditions the liquid is said to be moving as a ' forced vortex " with velocity q = co r at any radius 
r from the axis. [This contrasts with "free-vorte x" motion in which the fluid velocity varies 



With co constant, all fluid experiences an acceleration co 2 r (centrifugal) directed towards the axis 
of rotation, and for equilibrium of a typical small horizontal element: 


-z~or oA = p 8A Sr oi 2 r ... . (4.13) 
p a ) 2 r = p 


In the vertical direction, the usual expression for pressure distribution in a static fluid holds: 



Clearly, in such circumstances, the pressure is varying with both r and y: p = f(r, y),and 


dp = dr + .... (4.16) 

dr dv v y 


Hence: dp = p orr dr - pg dy, and if using a liquid with constant □, this can be integrated: 


p = - p co 2 r 2 — pgy + constant ... . (4.17) 
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• This reveals that all isobars in such a rotating liquid are paraboloids with the form: r 2 = ky + constant 

• and the free surface of the liquid (being an isobar) will also take this form. 

• We can find the equation of each isobar lines by applying the boundary condition of the center of 
rotation which have the coordinate of (0, y 0 ). 


<£xanwle (4) : The coffee up in example 1 is removed from the drag racer, placed on a turntable, 
and rotated about its central axis until a rigt ~ ^ ~V mode occurs. Find (a) the angular velocity 

which will cause the coffee to just reach the\ I the cun and (b) the gage pressure at point A 

for this condition. 


Solution : 



The cup contains 7 cm of coffee. The remaining distance of 3 cm up to the lip must equal the 
distance h/2 in Fig. 2.23. Thus 


h n _ il-R- a 2 (0.03 m) 2 
2 " 0 03 m " ~4f- ~ 4(9.81 m/s 2 ) 


Solving, we obtain 


ft 3 = 1 308 


or 


il = 36.2 rad/s = 345 r/min 


Arts, (a) 


To compute the pressure, it is convenient to put the origin of coordinates / and z at the bottom 
of the free-surface depression, as shown in Fig. E2.13. The gage pressure here is p 0 — 0, and 
point A is at (r, = (3 cm, ”4 cm). Equation (2.63) can then be evaluated 

p A = 0 - (1010 kg/m 3 )(9.81 m/s 3 )(-0.04 in) 

+ j( 1 0 1 0 kg/m 3 )(0.03 m) 2 (1 308 rad 2 /s 2 ) 

= 396 N/ in 2 + 594 N/m 2 = 990 Pa Arts, (b) 


This is about 43 percent greater than the still-water pressure p A = 694 Pa. 


example (5) : This tank is fitted with three piezometer columns of the same diameter. Calculate 
the pressure heads at points A & B when the tank is rotating at a constant speed of 100 rpm. 


Rotating stopped 


V 


I 


Pa — Vb 


100 rpm 
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Solution : 

1 58 1 

During rotation the surface in the piezometer coll] pust be on a parabola of constant pressure, 
and liquid lost from the central column must appear in the other two. Accordingly : 

PB-PA = (*0 2 W 2 = 171ft 
Y 2g 

^ = 5 - - (1.71) = 3.86 ft 

y 3 v y 

— = 5 + - (1.71) = 5.57 ft 

Y 3 v y 


@xattwle (6) : In a laboratory experiment, a circular tank of radius r 0 = 2 m is filled with water to 
an average depth d = 75 mm and rotated with an angular speed at = 5 rpm about its center. 
Find the shape of the free surface of the water when it is in rigid body rotation. 

Solution : 

The free surface is a parabola 
z - z 0 = co 2 r 2 / 2g 

with z 0 = the presently unknown depth at r = 0. The mass of fluid in the tank must be conserved 


so : 


pnr 0 2 d = p f* 0 2nr z dr = f^ 0 r 3 dr + 2npz 0 / Q to r dr 


/ 4 

2 i r o , 2 

r ° d = ~7i ~ + r ° z ° 

And : 


z 0 = d- = 0.075 - (— ) 2 — — = 0.047 m = 47 mm 

° 4 g v 60 J 4*9.81 

Thus : 

z = 0.047 + 0.014 r 2 < 0.103 = 103 mm 

that is, the water depth z ranges from 47 mm at the tank center to 103 mm at the tank wall. 
Compare these differences to those in the confined tank results above. Are the same principles 
operative? 
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Chapter Jftbe 

TYPES OF FLOW (FLUID FLOW) 


S. I. Introduction : 

Motion (flowing) of a fluid subjected to unbalanced forces or stresses that caused if the fluid mass 
was subjected to hydraulic gradient (e.g. tilting of free surface by certain angle or connect two 
containers have different levels). The motion continues as long as unbalanced forces are applied. 


\h 2 


FPi = pg./i2. 8 A 5 a 


hi 


Free surface 

FPi= pg.hi. 5A 

Flow 


FFfFF 2 

5.2. Ideal Fluid Flow J fUUUll UJ T7TTWC1 J 1 ^ 1 /UllllUll U UOUllFS E(]Ull tion)’ 

Consider a small element of fluid aligned along a streamline. It has a cross sectional area SA , 
pressure is assumed uniform across its ends SA, and the local velocity is defined v. 
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60 


Applying Newton's laws of motion to the flow through the element along the streamline: 


Theforce (in the direction ofmotion along the streamline) = mass x acceleration. 


1 . First we look at the acceleration of the fluid element. 

dv 

• Ignoring the possibility that the flow might be steady, — ^ 0 

dt 

• v can change with time t , and also with position s along the streamline. 

i.e. v=f(t,s). 


• Hence, if the element moves a distance ds in time St, then the total change in velocity 
Sv is given by: 



and in the limit as St tends to zero, the "substantive" derivative is given as: 


spatially temporarily 

r-H rS 


dv 

Sv 


dv , 

dv 



lim — = 

dv ds ^ dv 

v — + 


• • (5-2) 

dt 

8t — >0 

ds dt dt 

ds 

dt 


• for a steady flow the local velocity at a point ( 


dv 

~dt 


term under such circumstances will be zero. 


does not vary with time, so the last 


2. Looking now at the forces acting on the element and applying Newton's 2 nd law: 


p — ^ 5s) SA — pSA + pSA Ss g sin 6 = p Ss SAv ^ ... . (5.3) 


dividing through by SA. Ss and defining -Sz = Ss sin 0 , we have that: 



and in the limit as Ss tends to zero, 
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Or 


Euler’s equation (for ideal, steady flow) 


— + v dv + g dz = 0 ... . (5.7) 


This is a form of Euler's equation, and j 61 j p, v, and z in flow field. 

| I 

- it then becomes possible to integrate it - giving: 


p + - pv 2 + pgz = C .... (5.9) 
I pr 2,g v u 


Bernoulli’s equation (for ideal, steady flow) 


The three equations above are valid for incompressible, frictionless steady flow, and 
what they state is that total energy is conserved along a streamline. 

The first of these forms of the Bernoulli equation is a measure of energy per unit mass, 
the second of energy per unit volume, and the third of "head", equivalent to energy per unit 
weight. 

In the second equation, the term p is the static pressure, {Vipv 2 } is the dynamic pressure, 
pgz is the elevational term, and the SUM of all three is known as the stagnation (or total) 
pressure, po 

In the third equation, p/pg is known as the pressure head (or flow work head or flow energy 
head), which is the work done to move fluid against pressure, v 2 /2g as the kinetic head 
(dynamic energy head or velocity head), z is the potential head (elevation head) and the 
sum of the three terms as the Total Head H. The sum of first and third tern of 3 rd equation 
is called the piezometric head respect to piezometer's tube. 

where C is a constant along a streamline. 

> For the special case of irrotational flow, the constant C is the same everywhere in the flow 
field. 

> Therefore, the Bernoulli equation can be applied between any two points in the flow field if 
the flow is ideal, steady, incompressible, and irrotational. 
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> i.e. for two points 1 and 2 in the flow field: 


— + — + z-l = ^- + ^ + z 2 ... . (5.U) 

pg 2 g pg 2 g 


€xanwle (l): A n open tank filled with water at\ 62 ins through a port at the bottom of the 
tank. The elevation of the water in the tank Ts~IO m above the drain. The drain port is at 
atmospheric pressure. Find the velocity of the liquid in the drain port. 

Assumptions: 

1. Flow is steady. 

2. Viscous effects are unimportant. 

3. Velocity at liquid surface is much less than velocity in drain port. 



© 




i T 

1 
1 
1 
1 
1 
1 
\ 


&2g 

z \ 

\ 

\ 

\ 

\ 

\ 

\ 

< 

CD 


V 


: ► i 

I 

: ZJj 

r -2 

Streamline - 





Solution: 


1. The Bernoulli equation between points 1 and 2 on streamline: 


Pi | Vi 
P3 2 g 


Pi . iV 


+ Z 1 — h 


pg 2 g 


+ Z2 


2 . 


3 . 


The pressure at the outlet and the tank surface are the same (atmospheric), so pi = pi. The 
velocity at the tank surface is much less than in the drain port so Vi » Vi. Solution for 


Vi: 


0 + 0 + (zi — z 2 ) — 0 + — 

2 g 


Velocity calculation 

V 2 = j2g( z,-z 2 ) = V2 gH = a /2 x 9 . 8 1 x 10 = 14 m/s 
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Which states that the speed of efflux is equal to the speed of free fall from the surface of the 
reservoir. This is known as Torricelli’s theorem. 


! 63 | 

€xanwle (2): Calculate the flow rate through thtsyipeline and nozzle. Calculate the pressure at 
points 1, 2, 3 and 4 in the pipe and the elevation of the top of the jet’s trajectory. 



Solution: 

First sketch the energy line; it is evident that, at all points in the reservoir where the velocity is 
negligible, (z + p/y ) will be the same. Thus the energy line has the same elevation as the water 
surface. Next sketch the hydraulic grade line; this is coincident with the energy line in the 
reservoir where the velocity is negligible but drops below the energy line over the pipe entrance 
where the velocity is gained. The velocity in the 12- in. Pipe is everywhere the same, so the 
hydraulic grade line must be horizontal until the flow encounters the constriction upstream from 
section 2. Here, as velocity increases, the hydraulic grade line must fall (possibly to a level below 
the constriction). Downstream from the constriction the hydraulic grade line must rise to the 
original level over the 12- in. Pipe and continue at this level to a point over the base of the nozzle 
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at section 4. Over the nozzle the hydraulic grade line must fall to the nozzle tip and after that 
follow the jet, because the pressure in the jet is everywhere zero. 


Since the vertical distance between the energy line and the hydraulic grade line at any point is the 
velocity head at that point, it is evident that: 


100 

2g 

Thus : 


and therefore v 5 = 80*2. fvs { 

! 64 ; 

1 I 


2 

Q = 80.2 = 10.95 cfs 

From continuity considerations, 

"i = © 2 
And thus : 


Vl 


2 


2 9 



Similarly : 


Vjf_ 
2 9 



(^) 4 100 = 3/t 


15.3 ft 


Since the pressure heads are conventionally taken as the vertical distances between the pipe center 
line and hydraulic grade line, the pressure in pipe and constriction may be computed as follows: 


pi/y = 60 - 3 = 57 ft ; pi = 24.7 psi 

p 2 /y = 10 - 15.3 = -5.3 ft ; p 2 = 4.7 in. Hg vacuum 


p 3 /y = 40 - 3 = 37 ft ; p 3 = 16 psi 
pVy = 102 - 3 = 99 ft ; p 4 = 42.8 psi 

The velocity of the jet at the top of its trajectory (where there is no vertical component of velocity 
) is given by : 
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V6 = 80.2 cos 30° = 69.5 fps 
And the elevation here is : 

300 - (69.5) 2 / 2g = 225 ft. „ , 

! 65 ; 

€xamvle(3): The barometric pressure is 14 psia.- irhat diameter of constriction can be expected 
to produce incipient cavitation at the throat of the constriction ? 


E.L. 



Solution : 

From tables (appendix 2 : physical properties of water (Elementary fluid mechanics by Vennard 
and Streeter) 

y = 62 lb/ft 3 and p v / y = 2.2 ft ( at temperature of 100° F) 
p B / y = (14*144)/62 = 32.5 ft. 

Construct energy line and hydraulic grade line as indicated on the sketch. Evidently the velocity 
head at the constriction is given by : 

^- = 35 - 10 + 32.5 - 2.2 = 55.3 ft. 

2 9 1 

^ = 35 ft 

Since, from continuity considerations : 
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^6 


f. (jr=(f) 2 = 

V c o V c 


35 

55.3 


d c = 5.35 inch 


66 


5.3. Flow Classification (Flow patter/^ 1 


5.3.1. Uniform or Non Uniform Flow: 


Uniform flow is a flow in which the velocity does not change along a streamline, i.e. 



ds 


In uniform flows the streamlines are straight and parallel . 

Nonuniform (low is a flow in which the velocity changes along a streamline, i.e. 



ds 


5.3.2. Steady or Unsteady Flow: 

Steady flow: is a flow in which the velocity at a given point on a streamline does not change with 

dv _ 

time: — = 0 

dt 

dv 

Unsteady flow exists if: — ^ 0 

dt 

Combining the above we can classify any flow in to one of four types: 

• Steady uniform flow. Conditions do not change with position in the stream or with time. An 
example is the flow of water in a pipe of constant diameter at constant velocity. 

• Steady non-uniform flow. Conditions change from point to point in the stream but do not 
change with time. An example is flow in a tapering pipe with constant velocity at the inlet - 
velocity will change as you move along the length of the pipe toward the exit. 

• Unsteady uniform flow. At a given instant in time the conditions at every point are the same, 
but will change with time. An example is a pipe of constant diameter connected to a pump 
pumping at a constant rate which is then switched off. 

• Unsteady non-uniform flow. Every condition of the flow may change from point to point 
and with time at every point. An example is surface waves in an open channel. 
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5.4. Control Volume > 67 ! 


A control volume is a finite region, chosen by the analyst for a particular problem, with open 
boundaries through which mass, momentum, and energy are allowed to cross. The analyst makes 
a budget, or balance, between the incoming and 


In 

1 

1 

l 

1 

■1 

Control 

1 

1 

! 


1 

1 

1 

1 

Volume 

i m 

i 

i 

i 

i 


outgoing fluid and the resultant changes within the control volume. Therefore one can calculate 
the gross properties (net force, total power output, total heat transfer, etc.) with this method. 

With this method, however, we do not care about the details inside the control volume (In other 
words we can treat the control volume as a "black box.") 

5.5. Flow Kale ( Discharge) 


Mass Flow Rate = — : : = p . volume of flow rate = Discharge ( Q ) ... . (5.12) 

Time, taken to accumulative this mas£ __ _ 1 


Volume flow rate (Discharge) : 


Volume Flow Rate = = 

Time taken to accumiiMixe this, mlume 1 


If we need to calculate the discharge by knowing the cross sectional area of a pipe is A and 
Also we must know the Mean velocity is V m . 


Q = A . V m .... (5.14) 
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5.6. Continuity Equation (Conservation of Mass ) 


Mass entering per unit time = Mass leaving per unit time + Change of mass in control volume per unit time 



Applying to a stream tube as shown in figure below : 



Mass entering per unit time = Mass leaving per unit time 


piSAi Vi = p 2 5A 2 V 2 .... (5.15) 


or for steady flow : 



Continuity equation 
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to explain all above principles , the following example will show how to compute the inflow 
and outflow 


1. What is the outflow? 


control volume 



1.5 m /s 


Qir» — Qout 


1.5 + 1.5 = 3 
G out = 3.0 m 3 /s 


2. What is the inflow? 


u = 1 .5 m/s 
A = 0.5 m 2 


u = 1.0 m/s 
A = 0.7 m 2 


i r 


Q = 2 8 m /s 


] 


u = 0.2 m/s 

1 


A = 1.3 m 2 




h 1 

Q 


Q - Area x Mean Velocity = Au 

Q + 1 .5*0.5 + 1*0.7 = 0.2x1. 3 + 2.8 
Q = 3.72 m3/s 


Example 14): Calculate the flow rate of gasoline (s.g. =0.82) through this pipe line, using first 
the gage readings and then the manometer reading ? 

Solution : 


With sizes of pipe and constriction given, the problem centers around the computation of the 
value of the quantity (pi/ y + zi - pi/ y - zi ) of the Bernoulli equation. 


Taking datum at the lower gage and using the gage reading, 




20*144 10*144 

0.82*62.4 0.82*62.4 


4 = 24.2 ft of gasoline 


To use the manometer reading, construct the hydraulic grade line levels at 1 and 2. It is evident 
at once that the difference in these levels is the quantity (pi/ y + zi -pil y-zi) and, visualizing 
(pi/ y ) and (p 2 / y) as liquid columns, it is also apparent that the difference in the hydraulic grade 
line levels is equivalent to the manometer reading. Therefore : 
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fPi , V 2 \ 18.7 r 13.57—0.82^ - f r 

(j + z 1 -j-z 2 ) = —( — ) = 24.2 ft of gasoline 

Insertion of 24.2 ft into the Bernoulli equation -followed by its simultaneous solution with the 
continuity equation yields a flowrate Q of 8.(J 70 ; 


E.L. 



<£xmtwle (5) : (Venturi Section} Piezometric tubes are tapped into a Venturi section as shown in 
the figure. The liquid is incompressible. The upstream piezometric head is 1 m, and the 
piezometric head at the throat is 0.5 m. The throat section is twice large as in the approach 
section. Find the velocity in the throat section. 
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Problem Definition 



Situation: Incompressible flow in venturi section 71 ^metric heads and sections ratio given. 
Find: Velocity (in m/s) in venturi section. 

Assumptions: Viscosity effects are negligible, and the Bernoulli equation is applicable. 


Plan: 


1. Write out the continuity equation, incorporating sections ratio: 

2. Write out the Bernoulli equation, incorporating continuity eq. and solve for throat velocity. 

3. Substitute in piezometric heads to calculate throat velocity. 

Solution : 

The continuity equation Bernoulli gives: 


vA = vA 

A= 2A 

v 2 = 2v, 

The Bernoulli equation with V 2 = 2vi gives 


Pg 2g 1 pg 2 g 
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-^ + Zl -(-^ + z 2 ) = 

pg pg 


v„ 


V, 


2g 2g 


1. 0-0.5 = 


xv, 


v i = 


2g 


2g 


x0.5 =1.81 m/s 


v„ =2x1.81 = 3.62 m/s 


2 g 


5. 7. Applications of Bernoulli’s and Cp^^uity Equations : 

5.7. J. Flow_ fro/// Ta/iAs (F/ow T/iroite/i A Smalf_ Orifice) 

Flow from a tank through a hole in the side : 




A 


h 


V 



The edges of the hole are sharp to minimize frictional losses by minimizing the contact between 
the hole and the liquid. 

The streamlines at the orifice contract reducing the area of flow. 

This contraction is called the vena contracta 


The amount of contraction must be known to calculate the flow 
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Apply Bernoulli’s eq. along the streamline joining point 1 on the surface to point 2 at the centre 
of the orifice. 

At the surface velocity is negligible (ui = 0) and the pressure atmospheric (pi = 0). 

At the orifice the jet is open to the air so again the pressure is atmospheric (p 2 = 0). 

If we take the datum line through the orifice, then zi = h and Z 2 =0, leaving 


h = 


v 2 

2 ~9 ' 


and v 2 = y[2gh 


This theoretical value of velocity is an overestimate as friction losses have not been taken into 
account. 

A coefficient of velocity is used to correct the theoretical velocity, 

V actual = Cv V theoretical 

* 73 i 

Each orifice has its own coefficient of veloci] y usually lie in the range (0.97 - 0.99) 

The discharge through the orifice = jet area x jet velocity 

The area of the jet is the area of the vena contracta not the area of the orifice. 

We use a coefficient of contraction to get the area of the jet 

A actual = Co A orifice 

Giving discharge through the orifice: 

Q = vA 

Qactual = A actual V actual 

= Cc Cv A orifice V theoretical 


— Cd Aorifice V theoretical 


Qactual — Cd Aorifice <\j 2 .... (5.17) 


Where Cd is the coefficient of discharge , 


C d = C c * C v .... (5.18) 
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5.7.2. Time_ wanted to ex/iaiist t/ie water ofta/ti : 

We have an expression for the discharge from the tank 


Q= C d Ao J2gh ... . (5.18) 





St = ^ ... . (5.21) 

C d A 0 J 2 g Vh v ' 


hi 


We can use this to calculate how long it will take for level in the to fall 
As the tank empties the level of water falls. 

The discharge will also drop. 

The tank has a cross sectional area of A. 

In a time St the level falls by Sh 
The flow out of the tank is : 

p j j ' " Q (-ve sign as Sh\ ling) 

This Q is the same as the flow out of the orifice so : 




Integrating between the initial level, hi, and final level, In, gives the time it takes to fall this height 


[ CgAJTl^^- V^] - ■ (5.22) 


5.7.3. P/tot Static Tube ; 

Two piezometers, one as normal and one as a Pitot tube within the pipe can be used as shown 
below to measure velocity of flow. 
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By applying Bernoulli’s eq., we have the equation for pi : 

1 1 

Pi = pi + ~P Vi 2 i=> pgh2 - pghi + -pvi 2 



We now have an expression for velocity from two pressure measurements and the application of 
the Bernoulli equation. 

The necessity of two piezometers makes this arrangement uncomfortable. 

The Pitot static tube combines the tubes and they can then be easily connected to a manometer. 



The holes on the side connect to one side of a manometer, while the central hole connects to the 
other side of the manometer 

Using the theory of the manometer, 

Pa= pi+ pg (X-h) + pmangh and p B =p 2 +pgX 
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PA = Pb pi+ pg (X-h) + pmangh =p 2 +pgX 

we know that p 2 = pi + 0.5 pvi 2 , giving : 

pi + hg (pman-p) = pi + 0.5p Vl 2 



The Pitot/Pitot-static is: 

• Simple to use (and analyse) 

• Gives velocities (not discharge) 

• May block easily as the holes are small. 


5.7.4. Ve/ltl/n Meter : The Venturi meter is a device for measuring discharge in a pipe. It is 
a rapidly converging section which increases the velocity of flow and hence reduces the 
pressure. It then returns to the origii" - ^ "Imensions of the pipe by a gently diverging 
‘diffuser’ section. ! ! 



Apply Bernoulli along the streamline from point 1 to point 2 (eq. 5.1 1): 


Ei + V = Ei + v -l + Z2 

pg 2 g pg 2g 


By continuity equation : 
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Q = Vi Ai = V2 A2 '=> V>2 — —— 
A 2 


Substituting and rearranging gives : 



Note The theoretical (ideal) discharge is v* A and actual discharge takes into account the losses 
due to friction. We include a coefficient of discharge (Cd ~ 0.9) 


Q ideal = VlAl 

Qactual— Cd Qideal — Cd VlAl 


Qactual — CdA\A2 


2g[ ElZP2 + 

L pg 


In terms of the manometer readings : 
pi+pgzi =p 2 +p man gh +pg(z 2 - h) 
Pl ~ P2 + Zl - z 2 = h(^-~ 1) 



2gh[ Pman- !] 

^l 2 -^2 2 
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This expression does not include any elevation terms, (zi or zi) 

When used with a manometer, the Venturimeter can be used without knowing its angle. 


Ve/ztzzrimeter design^ 

• The diffuser assures a gradual and steady deceleration after the throat. So that pressure 
rises to something near that before the meter. 

• The angle of the diffuser is usually between 6 and 8 degrees. 
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• Wider and the flow might separate from the walls increasing energy loss. 

• If the angle is less the meter becomes very long and pressure losses again become 
significant. 

• The efficiency of the diffuser of increasing pressure back to the original is rarely greater 


• Care must be taken when connecting the manometer so that no burrs are present. 

S. 7. 5. Notches qnd_ Weirs_ : 

• A notch is an opening in the side of a tank or reservoir. 

• It is a device for measuring discharge 

• A weir is a notch on a larger scale - usually found in rivers. 

• It is used as both a discharge measuring device and a device to raise water levels. 

• There are many different designs of weir. 

• We will look at sharp crested weirs. 

5.7.5.1. We/r Assumptions ; i"”"! 

• velocity of the fluid approaching the weir is small so we can ignore kinetic energy. 

• The velocity in the flow depends only on the depth below the free surface v = yflgh. 
These assumptions are fine for tanks with notches or reservoirs with weirs, in rivers with high 
velocity approaching the weir is substantial the kinetic energy must be taken into account. 

5.7.5.2. General_ Weir Equatio/i : 

Consider a horizontal strip of width b, depth h below the free surface 


than 80%. 



velocity through the strip, v = ^2gh and discharge through the strip, 
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8Q = Av = b 8h ^Jlgh 

Integrating from the free surface, h = 0, to the weir crest . H = H, gives the total theoretical 
discharge 


H 1 / 

Qtheoretical ~ /q 2 dh ■■■ ■ (5-29) 


This is different for every differently shaped weir or notch. We need an expression relating the 
width of flow across the weir to the depth below the free surface 

5. 7.5.3. Recta/igu/ar Weir : 

The width does not change with depth so ( b = constant = B ) and substituting this into the general 
weir equation gives : 


Qtheoretical = B /" 6ftV 2 dh = 2 -B ... . (5.30) 


To get the actual discharge we introduce a coefficient of discharge, Ca, to account for losses at 
the edges of the weir and contractions in the area of flow : 


Qactual = cABjlgH 3 / 2 ....(5.31) 


€xamvle(6) (Rectangular Weir} : Water enters the. Mijlwood flood storage area via a rectangular 
weir when the river height exceeds the weir cj 79 wr design purposes a flow rate of 162 liters 
/ s over the weir can be assumed 

1. Assuming a height over the crest of 20cm and Cd=0.2, what is the necessary width , B, of 
the weir? 

2. What will be the velocity over the weir at this design? 


57.5.3. ‘V’ Notch Weir: 

The relationship between width and depth is dependent on the angle 0 of the “V”. 



The width, b, a depth h from the free surface is : 
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B = 2 (H-h)tan(^) 

So the discharge is : 

Qtheoretical = 2 j2g tan (|) /"(// - h)h 1/2 dh = 2^2 g tan (0 Hh 3/2 - ^ h 5/2 ]% 


Qtheoretical = ^ ^ tan (f) H S / 2 ... . (5.32) 


The actual discharge is obtained by introducing a coefficient of discharge : 


Q actual = Cd tan d)tf 5/2 -. (5.33) 


€xample(7) ‘V’ Notch Weir Water is flowing over a 90o ‘V’ Notch weir into a tank with a cross- 
sectional area of 0.6m2. After 30s the depth of the water in the tank is 1.5m. If the discharge 
coefficient for the weir is 0.8, what is the height of the water above the weir? 


5. S. The Momentum Equation : \ ! 

is a statement of Newton’s Second Law. It relates the sum of the forces to the acceleration or 
rate of change of momentum. From solid mechanics you will recognize 


What mass of moving fluid we should use? 

We use a different form of the equation. Consider a stream tube and assume steady non-uniform 
flow: 
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In time St a volume of the fluid moves from the inlet a distance uiSt, so 

volume entering the stream tube = area x distance =AiviSt 
mass entering stream tube = volume density =pAivi St 
momentum entering stream tube = mass velocity = p A1V1 Stvi 
Similarly, at the exit, we get the expression: 
momentum leaving stream tube = p A2U2 St 112 
By another reading of Newton’s 2 nd Law. 
where Momentum = m * v 


Force = mass x acceleration 


dv dmv 


rate of change of momentum 


dt dt 


F = 


lp2^2 v 2 St v 2 ~ P±A^V ifitVi) 
St 


We know from continuity that: Q = Ai vi=A2V2 
And if we have a fluid of constant density, 


81 ! 



The Momentum equation 
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This force acts on the fluid in the direction of the flow of the fluid 





The previous analysis assumed the inlet and outlet velocities in the same direction (i.e. a one 
dimensional system). 

What happens when this is not the case? 

We consider the forces by resolving in the directions of the co-ordinate axes. The force in the x- 
direction : 

Fx = pQ (V2 cos 62 - Vl COS 61) = pQ (V2x - Vlx ) 

And the force in the y-direction: 

F y = pQ (v2 sin 62 - vi sin 61) = pQ (vi y - vi y ) 

The resultant force can be found by combining these components 


Resultant 


F X 2 + Fy 2 


(5.35) 


And the angle of this force : 


6 = tan 1 (— ) ... . (5.36) 

SixZ 



This hydrodynamic force is made up of three components: 

Fr = Force exerted on the fluid by any solid body touching the control volume 

Fb = Force exerted on the fluid body (e.g. gravity) 

I 82 i 

Fp = Force exerted on the fluid by fluid pre^m-v-butside the control volume 
So we say that the total force, Ft, is given by the sum of these forces: 
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Ft = Fr + Fb + Fp 

The force exerted by the fluid on the solid body touching the control volume is opposite to Fr 
( action force). 

So the reaction force, R, is given by 

R = -Fr 


5. 8.1. Application of the Momentum Equation: 

5.8.1. 1. Forces on a Bend : 

Consider a converging or diverging pipe bend lying in the vertical or horizontal plane turning 
through an angle of 6. 

Here is a diagram of a diverging pipe bend 

Why do we want to know the forces here? 

As the fluid changes direction a force will act on the bend. 


This force can be very large in the case of water supply 
pipes. The bend must be held in place to prevent breakage at 
the joints. 

We need to know how much force a support (thrust block) must withstand. 
Step in Analysis: 

1 .Draw a control volume 

2. Decide on co-ordinate axis system 

3. Calculate the total force 

4. Calculate the pressure force 

5. Calculate the body force 

6. Calculate the resultant force. „ , 

! 83 ! 
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example 18): ( Forces on a Bend ) The outlet pipe from a pump is a bend of 45o rising in the 


100kN/m2 and the flow of water through the pipe is 0.3m3/s. The volume of the pipe is 0.07 5m3. 
[Ans. 13.95UN at 67 o 39* to the horizontal] 

Solution: 

1&2 Draw the control volume and the axis 
System 


vertical plane (i.e. and internal angle of 13 5o). The bend is 150mm diameter at its inlet and 
300mm diameter at its outlet. The pipe axis at the inlet is horizontal and at the outlet it is lm 
higher. By neglecting friction, calculate the force and its direction if the inlet pressure is 


Pi = 100 kN/m 2 
Q = 0.3 m 3 /s 
0= 45° 


A n = 0.177 m 2 As = 0.0707 m 2 


di = 0.15 m 


d 2 - 0.3 m 



3. Calculate the total force 


in the x direction 


F r x = 

- pQ(u 2 cos^-w,) 


by continuity 


i4,W| = A 2 u 2 = Q, so 


0.3 


1 6.98 m / s 


V\ = 



— - — = 4.24 m / .v 
0.0707 


I 


F Tx = lOOOx 0.3(4.24cos45- 16.98) 
= -41 93.68 N 


84 
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and in the y-direction 



- pQ{it 2 sin 0-0) 

= 1 000 x 03(4.24 s in 45) 
= 899.44 N 


4. Calculate the pressure force. 

F p = pressure force at ] - pressure force at 2 
&i = 0 , e-^e 

F p r - P\A\ cosO- p 2 A 2 co <s& = p i A l - p 2 Aj cos 8 
Fp, ~ i o,A, sinO- PjAj sin 0 = -p 2 A 2 sirtfl 


We know pressure at the inlet, but not at the outlet, we can use the Bernoulli equation to calculate 
this unknown pressure. 


Pi 

PZ 2 # pg 



f z 2 


The height of the pipe at the outlet is lm above the inlet. 

Taking the inlet level as the datum: 

zi = 0 , Z2 = lm , So the Bernoulli equation becomes: 


85 
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100000 16.98 2 +0= Pz + 4.24 s + ] 

1000 x 9.81 2x9.81 _ 1000x9.81 2x9.81 

p 2 =2253(51.4 N fm 1 

F Pj = 100000 x 0.01 77 - 22530 L4cos45 x 0.0707 
= 1770-11 266.34 = -949637 kN 


F Py = -225361.4 sin 45 x 0.0707 
= -11266.37 


5. Calculate the body force 

The body force is the force due to gravity. That 

is the weight acting in the -ve y direction. 

F b — —pg x volume 
= -1000x9.81x0.075 


F By = -735.75 N 


There are no body forces in the x direction, 


^.,=o 


6. Calculate the resultant force 


F = F 

1 r.v 


Ty 


^Ry + Py + ^ By 


Fr t — F Tx F Px F Bx 
= —41 93-6+ 9496.37 
= 5302.7 N 


F Sy — Fry ^Py F By 

= 899.44 + 1 1 266.3 7 + 735.7 5 
= 1 2901.56 vV 


86 ! 
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And the resultant force by the fluid is given by : 



= 1 3.95 kN 


And the direction of application is 



= 67.66° =67 *39’ 

The reaction force by the bend is the same magnitude but in the opposite direction 

R = -F it =-13.95 kN 


dxmtwle 19): The 6-cm-diameter 20° C water jet in Fig. strikes a plate containing a hole of 4-cm 
diameter. Part of the jet passes through the hole, and 

part is deflected. Determine the horizontal force n Plate 

required to hold the plat 


Solution : 


25 m/s 


25 m/s 
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Q = -( 0. 06 / 25 = 0. 0707 mfs 

4 

&„«, = ^0.04 /25 = 0.0314m 3 A 

F= PQ(K„-VJ 

for devided or branched flow 



F x = 99 8f 0 + 0 + 0.03 14 x 25 - 0.0707 x 25 ) 
F x = - 980 N 
R x = 980N 

r = o 


€xattwle (10): Water at 20° C exits to the standard sea-level atmosphere through the split nozzle 
in Fig. Duct areas are A1 - 0.02 m2 and A2 = A3 = 0.008 m2. If pi- 135 kPa (absolute) and 
the flow rate is Q2 - 03 = 275 m3/h, compute the force on the flange bolts at section 1. 
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Solution: With the known flow rates, wc can compute the various velocities: 

„ 275/3600 nv Vs n rr m „ 550/3600 m 

V, = V,= ^ — = 9.55 — ; V,= = 7.64 — 

0.008 nr s 0.02 s 

The CV encloses the split nozzle and cuts through the flange. The balance of forces is 

XF S =-W + Pi cas * A i =PQ,(-V 2 cos 30 o ) + / ?Q,(-V,cos30°)-/>Q i (+V 1 ), 




or: fU= 2(998) 


275 


(9.55cos30°) 


+ 998| - 


550 


V- 


3600 

= 1261 + 1 165 + 673 = 3100 N Am 


\3600 


(7.64) + (135000- 101350)(0.02) 


example (ll): A water jet of velocity V) impinges normal to a flat plate which moves to the right 
at velocity V c , as shown in Fig. Find the force required to keep the plate moving at constant 
velocity if the jet density is 1000 kg/m3, the jet area is 3 cm 2 , and Vj and Vc are 20 and 15 m/s, 
respectively. Neglect the weight of the jet and plate, and assume steady flow with respect to the 
moving plate with the jet splitting into an equal upward and downward half -jet. 


S' 

p = Po ! ■_ i cs 



Nozzle 


1/ 


n E 


Solution: 

For moving control volume with l/=l/cwe have 

V. =V -V =20-15 = 5 m/s 

in j c / 


CS 


v r v r 


A i(D 








By continuity equation we have: 
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Q =Q 


1 


AjV^W + Ay,, A l =A 2 = — A j 

V in = ^ V x + ^ V 2 , but from symmetry and neglacting the wight : V 1 — V 2 

v =v =v 

V in V 1 v 2 


F = pQ(V out -V in ) 

for devided or branched flow 

f = p(lfQ ,v , -Xg v ) 
f = p(y\A ,v y , -Ta v v ) 

F x = p( 0 — 1000 x 0.0003 x 5 x 5 j 
F x = —1 .5 N 

F x = F px + F l<x 

F Rx =F x - F px = —7.5 —0 — -7.5 
R X =-7.5N 
F — 0 


Example (12): The sluice gate in Fig. can control and measure flow in open channels. At sections 
1 and 2, the flow is uniform and the pressure is hydrostatic. The channel width is b into the 
paper. Neglecting bottom friction, derive an expression for the force F required to hold the 
gate. For what condition h2 is the force largest? For very low velocity V i 2 « ghl, for what 
value of h2/hl will the force be one-half of the maximum? 


1 


h | — - V, 


Sluice 

gate, width b 


r 


, A 7 
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Solution: The CV encloses the inlet and exit and the whole gate, as shown. From mass 

conservation, the velocities arc related by 

V l h 1 b = V 2 h 2 b, or: V^V^h^) 


The bottom pressures at sections I&.2 equal pg\\\ and jbghj, respectively. The horizontal 
force balance is 


TF k - -F gatc + ^ pgh, (h, b) ~ /?gh 2 (hib) = m(V 2 - V, ), m = /^bV, 


Solve for F pate =-pgbhai-(h,/h 1 )-]-/?h I bVi 


i-1 

h. 


Ans. 


For everything held constant except hi, the maximum force occurs when 


dF 

dlu 


= 0 


which yields 




yL 

§ h i 


\1/3 


Ans. 


Finally, for very low velocity, only the first term holds: F = ( l/2)/?gb (hf -hj). In this 
ease the maximum force occurs when hi - 0, or F maK - (l/2)pgbhf. (Clearly this is the 


special ease of the earlier results for F milx when V i =0.) Then for this latter ease of very 
low velocity, 

F = (l(2)F maK when h 2 =h,/V2. Ans. 


example ( 13 ) : Five liters/s of water enter the rotor shown in Fig. along the axis of rotation. The 
cross-sectional area of each of the three nozzle exits normal to the relative velocity is 18 
mm2. How large is the resisting torque required to hold the rotor stationary? How fast will 
the rotor spin steadily if the resisting torque is reduced to zero and (a) 0— 0 (b) 0— 30°, (c) 
0=60°? ; 


r = 0.h m 


Nozzle escit area normal to 
relative velocity - 16 rnm ? 
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example (MX : In the eductor shown in the following figure, water is injected through a 2 -cm 
nozzle at a speed of 30 m/s. The flow of water in the 5 -cm duct is 5 m/s. If the pressure 
downstream where the flow is totally mixed is 100 kPa, what is the pressure where the water is 
injected through the nozzle? Neglect the friction on the walls. 


p=100 kPa 


2cmJ 


-+-V=30 m/s 


*V=5 m/s 


5 cm 


Solution. 

Draw a control surface as shown with the appropriate force and momentum dia- 
grams. 



FD MD 


The velocity at the outlet may be obtained from the continuity equation for a steady 
flow. The mass How in the high-speed jet is 

ritjet = 1000 x ^ x 30 x 0.02 2 
= 9.42 kg/s 

The mass flow through the outer annular region is 

rn a = 1000 x jxSx (0.05 2 - 0.02 2 ) 

- 8.24 kg/s 

The velocity at the outlet is 

K — (™ jet + m a )/pA 

— (9.42 + 8.24)/(1000 x ^ x 0.05 2 ) 

= 9 m/s 
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The sum of the forces on the control surface is 

V F x = P\A - 1*2 A 

CS 

= ^ x 0.05 2 x (p! - 10 5 ) 

= 0.00106 x (pt - 10 5 ) N 

The control surface is not, accelerating and the flow is steady, so the momentum 
flux is 

" rhoV ] 'ritiVi^x = (0.12 +- 8.24) x 0 — 0.42 x 30 — 8.24 x 5 

cs OS 

— —164.86 kg - m/s * 

Equating the force and momentum flux 

0.00106 x (pi - 10 5 ) - -164.86 

= 15.0 kPa 
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Chapter &tx 

REAL FLOW IN CLOSED PIPES 


6.1. Introduction : 

The flow of a real fluid is vastly more complex than that of an ideal fluid, owing to phenomena 
caused by the existence of viscosity. Viscosity introduces resistance to motion by causing shear 
or friction forces between fluid particles and between these and boundary walls. For flow to take 
place, work must be done against these resistance forces, and in the process energy is converted 
into heat. 

6.2. Fundamental Equations : 


V 

h L 


V=Q/A 
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6.2.1. Reynolds Number : 

The effect of viscosity cause the flow of a real fluid to occur under three very different conditions 


1. Laminar flow 

2. Transitional flow 

3. Turbulent flow 


The characteristics of these regimes demonstrated by Reynolds number (R): 


^ = Reynolds number = ^ .... (6.1) 


(dimensionless value) 


In which : 

V = the mean velocity of the fluid in the pipe (m/s) 
d = diameter of the pipe (m) 
p = mass density of the fluid (kg/m 3 ) 
p = dynamic viscosity (N. s/m 1 2 or Pa .s or kg/m .s) 
v = kinematic viscosity (m 2 /s) 


Reynolds number for various flow regimes 

Flow Regime 

Reynolds Number 

Laminar 

<2000 

Transitional 

2000-4000 

Turbulent 

>4000 


1 . In laminar flow, agitation of fluid particles is of molecular nature only, and these particles 
are constrained to motion in essentially parallel paths by the action of viscosity. The 
shearing stress between adjacent moving layers is determined in laminar flow by the 

viscosity and is completely defined by the differential equation which was discussed in 

chapter one : 
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laminar flow occurs when Re < 2000. Laminar flow in a round tube is called Poiseuille flow or 
Hagen-Poiseuille flow in honor of pioneering researchers who studied low-speed flows in the 
1840s. For laminar flow only in pipes, shear -Stresses through flow is related by Newton law of 
viscosity ! 95 ; 



where y is the distance from the pipe wall. Change variables by letting y = ro - r, where ro is pipe 
radius and r is the radial coordinate. Next, use the chain rule of- calculus : 




dvy 


Rearrangement and substitutions with r =D/2 : 



dv = — r dr ... . (6.4) 

2uL ^ _ 


By integral with boundary condition of r = 0, V = V m ax 


^ ^max 


r 2 


2 u L 


(6.5) 


(parabolic velocity profile) 


From this equation we find that (at r=r 0 =D/2, V-0): 


I ^max 


h L y 


16/i L 


D 2 ....(6.6) 


The mean velocity^ V mean = Q /A .... (6.7) 


We find discharge by integration of V.dA overall the section area and we arrived to following 
result : 


V mean — V max — 


h L y 


32 uL 


D 2 .... (6.8) 


ft L = 32^r....(6.9) 


(Hagen-Poisseuille equation, 1839,40) 


By comparing this equation with Darcy equation, we find that: 
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For horizontal pipes: 


AP = (64/Re)(pV 2 /2)(L/D) = (64p/pVd)(pViMD) = 32pVL/d 2 .... (6.1 1) 


• Sometimes it’s more convenient to deal with volume flow rate (Q) rather than velocity (V). 
Thus we can write one last relation: 


AP = (1 28/7i) uQL/D 4 .... (6.12) 


(laminar flow only in horizontal pipes) 


Note the significance of this result: if you double the flow rate Q or the length of the pipe L, the 
pressure drop doubles (makes sense.) Also, for a given flow rate Q, if you double the diameter 
of the tube, the pressure drop decreases by a factor of 16! So use a little bigger pipe in your 
plumbing design. 


2. In turbulent flow, fluid particles do not remain in layers, but move in a heterogeneous 
fashion through the flow, sliding past other particles and colliding with some in an entirely 
haphazard manner that results in a rapid and continuous mixing of the fluid as flow occurs. 
These particles are observed to travel in randomly moving fluid masses of varying sizes 
called eddies. 



the friction factor depends not only on Re but also the roughness of the pipe wall, which is 
characterized by a roughness factor = e/d, where 8 is a measure of the roughness (i.e. height of 
the bumps on the wall) and d is (as always) the pipe diameter. The combined effects of roughness 
and Re are presented in terms of the Moody chart (R,e/D) still needs to be determined. For 
laminar flow, there is an exact solution for /since laminar pipe flow has an exact solution. For 
turbulent flow, approximate solution for/using log-law as per Moody diagram and discussed late 
or other approximations. 
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6.2.2. Darcy - Weisbaeh Equatio\ 97 ! 



In which : 

h l = head loss through a pipe due to friction effects (m) 

/= friction factor 

l = pipe length (m) 

g = gravitational acceleration (m / s 2 ) 

D - pipe diameter (m) 
v = fluid velocity (m/s) 

this equation, usually called the “Darcy equation” is still the basic equation for head loss 
caused by established pipe friction (not pipe fittings) in long, straight, uniform pipes. 

The equation of frictional stress is : 



To = F (v, d, p, p, e) 


And the velocity of friction may be expressed by : 



6.3. Equations wbicb relate between Reynolds number and frietion faetor (In ease 
of turbulent /low ) : 

/. 3/as/us Equation : 


0.316 

R 0.25 


R < 10 5 ....(6.16) 
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1. Burka Equation : 


f= S ... . (6.17) 


2. Von /Carman - Brandt/ Equations . 



iWIm 



J. Barr Equation . 


4. /Co ie brook- White Equation 


98 



njo 


5. Swamee Equation 


0.25 


e . 5.74' 17 ....(6.22) 


3 * 10 3 < R < 3 * 10 8 & 10" 6 < - < 2 * 10" 2 

~~ ~ D~ 



Reynolds number (Re) 


- Pipe Flow Problem Types 


Variable 

Type I 

Type II 

Type III 

a. Fluid 




Density 

Given 

Given 

Given 
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Viscosity 

Given 

Given 

Given 

b. Pipe 

Diameter 

Given 

Given 

Determine 

Length 

Given 

Given 

Given 

Roughness 

Given 

Given 

Given 

c. Flow 

Flowrate or Average 

Given 

Determine 

Given 

Velocity 





d. Pressure 

Pressure Drop or Determine Given Given 

Head Loss •- i 



x = turbulent shear stress 


s = eddy viscosity 


x ff = -JWxVy ■■■■(6.24) 


tr = Reynolds stress 

v x , v y = velocities of fluid particles in turbulent flow conditions at x and y directions respectively. 


v x v y = mean value of the product v x & v y 



In which : l = mixing length (mean distance in which the fluid particles transported by turbulence 
between layers) 

k = dimensionless constant of turbulence 


example (l): Water of kinematic viscosity 1.15*10' 6 m 2 /s flows in a cylindrical pipe of 30 mm 
diameter. Calculate the largest flow rate for which laminar flow can be expected. What is the 
equivalent flow rate for air ? 

Solution : taking R = 2100 as the conservative upper limit for laminar flow, 
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2100 = K ( °'° 3) ; V = 0.0805 mps 

Q water = 0.0805 * ^ (0.03) 2 = 5.69 * 10" 5 m 3 /s 

Vair = 1.37 * 10' 5 m 2 /s and Q a ir = 6.78*10' 4 m 3 /s □ ~ 12 Q wa ter 

exattwle (2): Water flows in a 150 mm diamei 100 \eline at a mean velocity of 4.5 mps. The head 
lost in 30 m of this pipe is measured experimentally and found to be 5.33 m. calculate the 
friction velocity. 

Solution : 

^0 4 S 2 

5.33 =f — — ; / = 0.026 

7 0.15 2 g J 

T , . - |o.026 „ _ , 

V 0 = 4.5 I— = 0.26 mps 


€xamvle (3): A turbulent flow of water occurs in a pipe of 2 ft diameter. The velocity profile is 
measured experimentally and found to be closely approximated by the equation v = 10 + 0.8 In 
y, in which v is in ft/sec and y ( the distance from the pipe wall )is in feet. The shearing stress 
in the fluid at a point 4 in. from the wall is calculated analytically from measurements of 
pressure drop to be 0.2 psf. Calculate the eddy viscosity, mixing length and turbulence constant 
at this point. 

Solution : 


dv 

dy 


— (10 + 0.8 In y) 


*5= ™=2.4 

y £) 


d 2 v _ d /0.8\ _ -0.8 _ -0.8 _ 
dy 2 dy \ y ) y 2 (Ay 

0.2 = s (2.4) ; s = 0.083 lb.sec / ft 2 
0.2= 1.94 / 2 (2.4) 2 ; / = 0.134 ft 


0.2= 1 - 9 ^ fc2( ^ 4)4 ; k = 0.401 

(— 7.2) 2 


0.2 = 1.94 k 2 (-) 2 (2.4) 2 


k = 


1= °- 401(2 - 4) = 0.134 ft 

(-7.2) 


0.401 
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1 = 0.401 (£) = 0.134 ft 

The magnitude of the eddy viscosity e when compared with the viscosity p (approximately 

0.000021b. sec / ft 2 ) is of special interest in that it provides a direct comparison between the large 
turbulent shear and small laminar shear fi'Yoi”' same velocity gradient. The mixing length, / , 

when compared with the pipe radius is foi ! be about 10 % of the latter dimension; this is a 

nominal value of correct order of magnitude, as is the turbulence constant, k . 


(Example (4) : (Head loss for laminar flow) Oil (S = 0.85) with a kinematic viscosity of 6 x 10-4 
mils flows in a 15 cm pipe at a rate of 0.020 m3/s. What is the head loss per 100 m length of 
pipe? 

Problem Definition 

1. Oil is flowing in a pipe at a flow rate of Q- 0.02 m3/s. 

2. Pipe diameter is D - 0.15 m. 

Find: Head loss (in meters ) for a pipe length of 100 m. 

Properties: Oil: S = 0.85, v = 6 x 10-4 m2/s. 

Solution : 


1 . Mean velocity 


y-2- 0.020 m 3 1 s _ 0.020 m^/s , 13mJ , s 

-4 " (r£5 a )M i((0 15m) a /4) 


2. Reynolds number 


E VD _ (1.13 a/ 0(0 15 m) _ =C3 

fix 10“*tn a f s 


3. Since Re < 2000, the flow is laminar. 


4. Head loss (laminar flow). 


. 3ZidV 32 tvLV _ S2i/L¥ 

f 1D 2 pgD 1 gD 3 

35(6* 10 ^ m ] t_ s)(100 m)qi3 m/s) 
(3.81 ml s 2 )(0. 15 m) s 


Review: Tip! An alternative way to calculate head loss for laminar flow is to use the Darcy - 
Weisbach equation as follows: 
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f 

k f 



_ 64 _ 

283 

2g 


9.S3 m 


0.226 

-0.226 


/ 1 00 m ' ! , 
I 0.15 m J 


I (1.13m /Q ^ 

^ 2 x 5 01 m ^ s 2 J 



gg a 


4 


=/iO)3tfJ UOJPUJ 


Moody Diagram 


Reynolds number R _ — consistent units 
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€xamvle (5) : ( Type 1 ) Determine the head lost to friction when water at 15 °C, v = 1.14 mm2 • 
s~l , flows through 300 m of 150 mm diameter galvanized steel pipe (s = 0.15 mm) at 50 L • s 


Solution : 

V = 50 x 10 _3 m 3 • s _1 /((7i/4)(0.15) 2 m 2 ) = 2.83 m- s _1 

Re = VD/v = 2.83 m- s” 1 x 0.15 m/(1.14 x 10" 6 m 2 • s” 1 ) = 3.72 x 10 5 

For galvanized steel, e /D = 0.001 

From Moody Diagram /= 0.0208, so 

H l = 4 x0.0208 x 300 m/(0.15 m) x (2.83 m- s" 1 ) 2 /(19.62 m- s" 2 ) = 67.89 m, 
say 68 m 


example 16): (Type II) Oil, with p =950 kg/m3 and v = 2 E-5 m2/s, flows through a 30- au- 
diometer pipe 100 m long with a head loss of 8 m. The roughness ratio is s/d = 0.0002. Find 
the average velocity and flow rate. 

Solution: 

By definition, the friction factor is known except for V: 


Guess / — 0.014, then V = V O. 4 7 1/0.014 = 5.80 m/s and Re t/ = Vd/v — 87,000. At Kc d = 
87,000 and dd = 0.0002, compute / new - 0.0195 [Eq. (6.64)]. 

New f~ 0.0195, V = V0.481/0.0195 = 4.91 m/s and Re t/ = Vd/v = 73,700. At Re^ = 
73,700 and dd = 0.0002, compute / ncw — 0.0201 [Eq. (6.64)]. 

Better/ = 0.0201, V = Vo.47 1/0.0201 = 4.84 m/s and Re, y = 72,600. At Re rf = 72,600 and 
dd = 0.0002, compute / ncw “ 0.0201 [Eq. (6.64)]. 

We have converged to three significant figures. Thus our iterative solution is 


1 



or /V 2 = 0.471 (SI units) 


To get started, we only need to guess/, compute V — V0 ,471/f, then get Re^, compute a better 
/ from the Moody chart, and repeat. The process converges fairly rapidly. A good first guess is 
the “fully rough” value for dd = 0.0002, or / — 0.014 from Fig. 6.13. The iteration would be as 


follows: 


V ~ 4.84 m/s 



Ans. 


104 


"D%. "padAil /iid dfl- 


"plaid "7KecA<zmc4 


Secmd Stacie 


€xamvle (7) : (Type III) Work previous Example backward , assuming that Q = 0.342 m3/s and s 
- 0.06 mm are known but that d (30 cm) is unknown. Recall L = 100 m, p = 950 kg/m3, v = 2 
E-5 m2/s, and hL = 8m 

Solution: 


First write the diameter in terms of the friction factor: 


tt (9.81 m/s 2 )(8 m)f/ 5 
_ 8 _ (100 m)(0.342 nrVs ) 2 


8 . 28 ^ 


or (/*» 0.655/ l/3 


( 1 ) 


in SI units. Also write the Reynolds number and roughness ratio in terms of the diameter: 


Re 4(0.342 nrVs) 21,800 

7t(2 E-5 m 2 /s )d d 

e _ 6 E-5 m 

d ~ d 

Guess/, compute d from (1), then compute Re^ from (2) and dd from (3), and compute a bet- 
ter/from the Moody chart or Eq. (6.64). Repeat until (fairly rapid) convergence. Having no ini- 
tial estimate for/, the writer guesses / = 0.03 (about in the middle of the turbulent portion of 
the Moody chart). The following calculations result: 


/« 0.03 d = 0.65 5(0.03) 1/5 = 0.325 m 


21,800 

Re., = 

» 67,000 

- = 1.85 E-4 

0.325 

d 

/ ncw - 0.0203 

then 

d ncvJ = 0.301 m 

Re^.ncw =72,500 

4 = 2.0 E-4 
a 

/better ** 0.0201 

and 

d = 0.300 m 


( 2 ) 

(3) 
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example 18) : SAE 30 oil at 20° C flows in the 3-cm-diameter pipe in Fig. which slopes at 37°. For 
the pressure measurements shown, determine (a) whether the flow is up or down and (b) the 
flow rate in m3/h 


Pb 


180000 


HGLr — * + z fj — 

pg 891(9.81) 


+ 15 = 35.6 m; HGL A =■ 


500000 


891(9.81) 

Since HGL a > HGL b the flow is up Ans.(a) 
The head loss is the difference between hydraulic grade levels: 


+ 0 = 57.2 m 


u 128//LQ 1 28(0.29)(25)Q 

h f = 57.2 - 35.6 = 2 1 .6 m = — -r- = r 

rrpgd 4 7T(89l)(9.8l)(0.03) 4 

Solve for Q - 0.0005 18 m Vs = 1.86 m^/h Ahs. (b) 
Finally, check Re - 4 pQ/( npd) *= 68 (OK, laminar flow). 


Example (0) : The velocity of oil (S = 0.8) through the 5 cm smooth pipe is 1.2 m/s. ere L = 12 m, 
zl = 1 m, z2 = 2 m, and the manometer deflection is 10 cm. Determine the flow direction, the 
resistance coefficient f, whether the flow is laminar or turbulent, and the viscosity of the oil. 



Solution: 

Based on the deflection on the manometer, the piezometric head on the right side of the pipe is 
larger than that on the left side. Since the velocity at 1 and 2 is the same, the energy at location 2 
is higher than the energy at location 1 . Since the a fluid will move from a location of high energy 
to a location of low energy, the flow is downward (from right to left). 

From energy principles: 
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+ z 7 + h. 
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Manometer equation 

pa -|- (2 m) Toil + (0 1 in) Toil ~ (01 in) Trig — t 1 111 ) Toil = Pi 

Algebra gives 


P 2 - P 1 
Toil 


P 2 Pj_ 

Toil 


- (2 in) — (0.1 ill) + (0.1 in) + ( lni ) 

Toil 

- (lm) + (01 in) - 1^ 

- C i m ) + (o.i m) 

0.6 m 


( 2 ) 


Substituting Eq. (2) into (1) gives 

(0.6 m) 


, , „ L V 2 

(-im) + /-— 

or 


/ = i 


<T)m 

- - m (^) 


/ 0.0908 


f 

64 

Ro 


0.0908 = 

64/t 


pVD 


or 


P - 

0 .0908/? V"Z> 
64 



0.0908 x (0.8 x 1000) x 1.2 x 0.05 


64 

= 

0.068 N-s/in 2 



Re 


VDp 

1.2 x 0.05 x (0.8 x 1000) 
0.008 


Thus, flow is laminar. 
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6.4. Minor Head Losses: 

For any pipe system, in addition to the Moody-type friction loss (referred as Major losses ) 
computed for the length of pipe, there are additional so-called minor losses due to: 

1. Pipe entrance or exit 

2. Sudden expansion or contraction 

3. Bends, elbows, tees, and other fittings 

4. Valves, open or partially closed 

5. Gradual expansions or contractions 

These losses may not be so minor in its effects; e.g., a partially closed valve can cause a greater 
pressure drop than a long pipe. 

Since the flow pattern in fittings and valves is quite complex, the theory is very weak. The 
losses are commonly measured experimentally and correlated with the pipe flow parameters. 

The data, especially for valves, are somewhat dependent upon the particular manufacturer’s 
design, so that the values listed here must be taken as average design estimates 

The measured minor loss is usually given as a ratio of the head loss hm through the device to 
the velocity head V 2 !{2g) of the associated piping system 


A single pipe system may have many minor losses. Since all are correlated with V 2 /( 2g), they 
can be summed into a single total system loss if the pipe has constant diameter 


Minor loss coefficient, K = — 

V / 2g 




K l = 2 based on 
pipe 

. velocity 


0. 50-in. 
diameter 


15 ft 


90° elbows 


l 
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L5 LOSS COEFFICIENTS FOR VARIOUS TRANSITIONS AND FITTINGS 


Description 

Sketch 

Additional Data 

K 

Source 

Pipe entrance 

\ V 


rid 


AT 

m r 

h L =KJ' Z /2g 

A 

4 I* 

0.0 


0.50 



/ “ 

T* 

0.1 


0.12 





>0.2 


0.03 


Contraction 

„ 

1 •? 


AT 

AT: 


h L = Kc^l / 2g 

* 


DiiD\ 

0 = 60° 

0 = 1 80° 

(10) 


i 

** 

0.00 

0.08 

0.50 



0.20 

0.40 

0.60 

0.80 

0.90 


0.08 

0.07 

0.06 

0.06 

0.06 


0.49 

0.42 

0.27 

0.20 

0.10 


Expansion 


h L = K^V2g 


th 




AT 

AT 

D\iD 2 

0 = 20° 

0 — 1 80^ 

0.00 


1.00 

0.20 

0.30 

0.87 

0.40 

0.25 

0.70 

0.60 

0 15 

0.41 

0.80 

0.10 

0.15 


90° miter bend 


90° smooth bend 



Without vanes Kb = 1.1 


With vanes 

AV 

rid 

1 

AT = 0.35 

2 

0.19 

4 

0.16 

6 

0.21 

8 

0.28 

10 

0.32 


(15) 


(15) 

(16) and (9) 


Th readed p i pe fill mgs G lobe v al ve — w ide ope n 
Angle val ve — wide open 
Gate valve — wide open 
Gate valve — half open 
Return bend 
Tee 

Straight-through flow 
Side-outlet flow 
90° elbow 


A v ^10.0 (15) 

AT = 5.0 
AT = 0.2 
A v = 5.6 
AT = 2.2 


A, - 0.4 
K,= 1.8 
AT = 0.9 
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45° elbow 


Kb = 0.4 
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€xamole (10) Find the distance through the pipeline in the figure below for H -10 m 
and determine the head loss Hfor Q - 60 L/s. 


Solution : the energy equation applied between points 1 and 2, including all the losses, 
can be written as : 


v 2 

H 1+U + 0=^-+U + 0+^+f 

2 g '» 


IV 2 


102 in V? 


V-? 


V? 


2 2g J 0.15 m 2 g 


+ 2(0.9)-^ + 10-^ 


2g 


2g 


in which the entrance loss coefficient is each elbow is 0.9, and the globe valve is 10. 
Then 


Wi = ^(13.3+680/) 

Zg 

When the head is given , this problem is solved as the second type of simple pipe problem. 
If / = 0.022, 


V 2 

10 = -^-[13.3 + 680(0.022)] 


and V '2 — 2.63 in/s. From Appendix C, 


1 / — 1.01 (J, m 2 /s 


= 0.0017 


(2.63 m/s) (0.15 m) 
_ 1.01 x 10~ e m 2 /s 


391,000 


H = 10 m 150-mm-diam clean cast-iron pipe £ 2 } 

* -a 



x 
T 

Standard elbows 
60 m — 


Globe valve 


Square- edged entrance 

Pipeline with minor losses. 


From Fig. 6.21, / = 0.023. Repeating the procedure gives V -2 = 2.60 in/s, R = 380,000, 
and / = 0.023. The discharge is 

Q = V 2 A 2 = (2.60 m/s) -^(0.15 m) 2 = 45.9 L/s 
For the second part, with Q known, the solution is straightforward: 

v * = %= Wi0^ = 3 - mm / s R = r>ur ’> uuo / = mixs 

and 


Hi 


(3.4 in/s) 2 
2(9.806 m/s 2 ) 


[13.3+680(0.023)] = 17.06 m 
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6.5. Pumps and Turbines: 

• Pumps and turbines represent the external source-sink term of energy head 

• Apply energy equation from the reservoir water surface to the outlet stream is: 

Zl 


In Ideal Flow: 


— + — + z 1 + hp = — + — + z 2 .... (6.26) 

y 2g 1 r y 2a z v y 


For a system with one size of pipe, this simplifies to : 




In real Flow with: 


Pl + ^ + Zl + hp = ^+ ^- + z 2 + £K^- + £^....( 6 . 28 ) 


V2 , V 2 


fL v 2 ‘ 


r 2 


y 2 


For a system with one size of pipe, this simplifies to : 


hp = (z 2 — Zl )+ ^(1+£K+ 2^). ...(6.29) 


Pump Power = pgQh .... (6.30) 


(hydraulic power of pump, Power out) 


Efficiency = -i^....(6.31) 

power j n power j n 
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example (11) : Water, p = 1000 kg/ m3 . and v = 1.14 X 10-6 m2/s, is pumped between two 
reservoirs at 6 l/s through 121 m of 50 mm-diameter pipe and several minor losses, as shown 
in Fig. The roughness ratio is s/d = 0.001. Compute the pump horsepower required. 



loss 


K 

Sharp entrant ) 


0.5 

Open globe valve (2 i 


6.9 

12-in ben 


0.15 

Regular 90° elbo 


0.95 

Half- do sal gate vah 

» 

2.7 

Sharp es 


1.0 

2 K= 12.2 
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6. 6. Pipes Systems : 




1 

t 


(ct> 


V i 





* 








f 1 *1 




U 2 (2) 



' i *i 




u 3 oo 

— ► ©3 



1 t 

1 


tw 


a. Pipes in Series: 


|hL A-B - hLl + hL2 + hL3 • ■ • • (6.33)| 


b. Pipes in Parallel: 




example 112): Given is a three-pipe series system , as in Fig. The total pressure drop is pA- pB - 
150,000 Pa, and the elevation drop is zA-zB = 5 m. The pipe data are in table below. The fluid 
is water, p = 1000 kg/m3 and v = 1.02 x 10 -6 m2/s. Calculate the flow rate Q in m3/h through 
the system. 


i 

A * 





•B 
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Pipe 

L, m 

di cm 

t/ mm 

dd 

1 

100 

8 

0.24 

0.003 

2 

150 

6 

0.12 

0.002 

3 

80 

4 

0.20 

0.005 


Solution: 


The total head loss across the system is 

a j _ Pa Pb I _ T = 150,000 

A ~* B pg <A B 1000(9.81) 

From the continuity relation (6.105) the velocities are 

v 2 =^v, =~-V x V^^Lv, = 4V, 

£"/ T 9 £l 3 


+ 5 m = 20.3 m 


and 


V^L 4 

Re^ - —2^ Re, - —Re, Re ? - 2 Re, 

V,r/, 3 


Neglecting minor losses and substituting into Eq, (6. 107), we obtain 


or 


A h A ^ B = 1^1250/, + ISOO^f Jf 2 + 2000(4) 2 / 3 J 
V? 

20.3 m = A-L ( ] 250 f, + 7900 f 2 + 32,000 f 3 ) 

-g 


( 1 ) 


This is the form which was hinted at in Eq. (6.108). It seems to be dominated by the third pipe 

loss 32,000/ 3 . Begin by estimating/1,/2, and / 3 from the Moody-ehart fully rough regime 

/, = 0.0262 f 2 = 0.0234 / 3 = 0.0304 

Substitute in Eq. (1) to find V 2 = 2#(20.3)/(33 + 185 + 973). The first estimate thus is V, = 
0.58 m/s, from which 

Re, = 45,400 Re 2 = 60,500 Re :? = 90,800 
Hence, from the Moody chart, 

/, = 0.0288 f 2 = 0.0260 / 3 = 0.0314 

Substitution into Eq. (1) gives the better estimate 

V, = 0.565 m/s Q = 4W7V, = 2.84 x 10“ 3 m 3 /s 

or Q, = 10.2 m 3 /h 

A second iteration gives Q = 10.22 m 3 /h, a negligible change. 


A ns. 
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example (13): Assume that the same three pipes in Example above are now in parallel with the 
same total head loss of 20.3 m. Compute the total flow rate Q, neglecting minor losses. 


Pipe 

L, m 

th c m 

mm 

dd 

1 

100 

8 

0.24 

0.003 

2 

150 

6 

0.12 

0.002 

3 

80 

4 

0.20 

0.005 


Solution: 


20,3 m = p 1250/, = p- 2500/. = p 2000/ 

-8 -8 


0 ) 


Guess fully rough flow in pipe 1:/, = 0.0262, V, - 3,49 m/s: hence Rej - V x d\tv - 273,000, 
From the Moody chart read/ = 0.0267; recompute V, = 3,46 m/s, (/, = 62,5 m 3 /h. [This prob- 
lem can also be solved from Eq, (6,66).| 

Next guess for pipe 2:f 2 ~ 0,0234, V 2 ~ 2.61 m/s: then Re 2 = 153,000, and hence / = 
0.0246, V 2 = 2.55 m/s, Q 2 = 25.9 m 3 /h. 

Finally guess for pipe 3:f~ 0,0304, V$ ~ 2,56 m/s; then Re 3 = 100,000, and hence / - 
0.0313, V 3 = 2.52 m/s, Q 2 = 1 1.4 m 3 /h. 

This is satisfactory convergence. The total flow rate is 

Q = Q\ + Qi + & = 625 + 25.9 + 1 1 .4 = 99.8 m 3 /h Am. 


ns ; 


i 


